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Abstract. We present exponential error estimates and demonstrate an algebraic convergence rate 
for the homogenization of level-set convex Hamilton- Jacobi equations in i.i.d. random environments, 
the first quantitative homogenization results for these equations in the stochastic setting. By taking 
advantage of a connection between the metric approach to homogenization and the theory of first- 
passage percolation, we obtain estimates on the fluctuations of the solutions to the approximate 
cell problem in the ballistic regime (away from flat spot of the effective Hamiltonian). In the sub- 
ballistic regime (on the flat spot), we show that the fluctuations are governed by an entirely different 
mechanism and the homogenization may proceed, without further assumptions, at an arbitrarily 
slow rate. We identify a necessary and sufficient condition on the law of the Hamiltonian for an 
algebraic rate of convergence to hold in the sub-ballistic regime and show, under this hypothesis, 
that the two rates may be merged to yield comprehensive error estimates and an algebraic rate of 
convergence for homogenization. 

Our methods are novel and quite different from the techniques employed in the periodic setting, 
although we benefit from previous works in both first-passage percolation and homogenization. The 
link between the rate of homogenization and the flat spot of the effective Hamiltonian, which is 
related to the nonexistence of correctors, is a purely random phenomenon observed here for the 
first time. 



1. Introduction 

We consider the Hamilton-Jacobi equation 

u\ + H(Du £ ,^,uj) = in R d x (0,oo), (1.1) 

where the Hamiltonian H = H(p,y,u) is level-set convex and coercive in p and depends on an 
element uj of an underlying probability space (f2,3~, P). If the action of translation on M rf is sta- 
tionary and ergodic with respect to the law of H, then, as e — > 0, the solutions u e = u £ (x,t,oj) 
of (jl.ip . subject to appropriate initial conditions, converge P-almost surely to the solution u of the 
deterministic equation 

u t + H(Du) = inR d x(0,oo) (1.2) 

with the same initial conditions, where the effective Hamiltonian H G C(M rf ) is level-set convex and 
coercive. This fundamental theorem concerning the qualitative theory of stochastic homogenization 
of Hamilton-Jacobi equations was proved for convex Hamiltonians by one of the authors [32] (see 
also Rezakhanlou and Tarver |31| ) and, more recently, by two of the authors [4] in the generality 
discussed here. 

In this paper, we present the first quantitative homogenization results for Hamilton-Jacobi 
equations in the stochastic setting. Throughout the paper we assume that the Hamiltonian 
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H = H(p,y,u) satisfies a finite range dependence hypothesis (a continuum analogue of "i.i.d.") in 
its spatial dependence. This essentially means that, for some fixed distance D > 0, the values of 
H(p, y, •) for y £ E are independent of those for y £ F provided that dist(£', F) > D. 

By a novel integration of probabilistic and pde techniques, we (i) obtain explicit estimates show- 
ing the probability that \u e (x,t,uj) — u(x,t)\ > A decays exponentially in A 2 , and (ii) identify a 
necessary and sufficient condition for the almost sure, local uniform convergence u £ — > u to proceed 
at an algebraic rate 0(s a ). The main results, including the precise assumptions, are stated in the 
next section. They essentially give the error estimates 



for the homogenization of (jl.ip . where the first inequality depends on a supplemental assumption 
on the law of H. 

The difficulty in obtaining estimates on the fluctuations of u £ (x, t, •) — u(x, t) is due in part to the 
fact that the dependence of u £ on H is highly singular. Understanding how the solutions depend 
on the random environment is very challenging. Difficulties of a similar nature occur, for example, 
in the theory of first passage percolation (see Kesten [18] and Alexander pQ) and in the study 
of the fluctuations of the Lyapunov exponents for Brownian motion in Poissonian potentials (see 
Sznitman [33] and Wiithrich [36J). As far as we know, the only previous result on the oscillations of 
solutions of Hamilton- Jacobi equations in random media is found in the work of Rezakhanlou [30J , 
who gave structural conditions on H in dimension d = 1 in which a central limit theorem holds. 
Such phenomena are not expected to appear in any dimension d>2 (see Remark 14. 2p . 

Our arguments rely crucially on adaptations of some of the probabilistic techniques of |18|. [I] , 
which are based on Azuma's inequality and the martingale method of bounded differences. This 
connection between first-passage percolation and the stochastic homogenization of (jl.ip . made 
explicit for the first time in this paper (as far as we know), arises naturally from an analogy 
between the passage time in percolation and solutions of the metric problem (see Remark IA.3I 
below). Using arguments inspired from |181 [I], we prove exponential error estimate and obtain 
rates of convergence for the homogenization of the metric problem. Then, by quantifying the new 
proof of homogenization recently introduced by two of the authors [1] , we transform the estimates 
for the metric problem into error estimates for the approximate cell problem (see (II. 4p below). 

The rate of convergence of periodic homogenization of Hamilton-Jacobi has been understood for 
some time and goes back to the work of Capuzzo-Dolcetta and Ishii [8], who proved that u £ and u 
differ by at most 0(e^). The periodic setting is much simpler to understand due to the fact that 
the cell problem has periodic solutions; that is, exact correctors exist. A quantitative version of 
the classical perturbed test function proof of homogenization due to Evans [12} [T3] then yields the 
convergence rate. Our main results stated in Section [2] do not encompass the periodic or almost 
periodic settings, since obviously an almost periodic function cannot be embedded into the random 
setting in such a way that it satisfies a finite range of dependence condition. However, as we show, 
our arguments yield a uniform rate of convergence in the almost periodic setting (see Section [7]) . 

In the stochastic environment, the situation is not only much more complicated but also qualita- 
tively different from the periodic setting. It is, therefore, necessary to devise a new strategy, since 
the usual proof of periodic homogenization, which is based on exact correctors and can be quantified 
to yield a rate, does not generalize to random enviroments. Indeed, as Lions and Souganidis [22] 
demonstrated with an explicit example, exact correctors do not exist, in general, for stochastic 
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Hamiltonians. The only known proofs of the qualitative homogenization of Hamilton-Jacobi equa- 
tions in the stationary ergodic setting are based on an application of the subadditive ergodic theo- 
rem to certain subadditive quantities (e.g., the m^'s below) and then showing that these quantities 
control, in an appropriate way, the solutions of (jl.ip . In order to obtain a convergence rate for 
the homogenization, one is therefore left with the twofold task of quantifying both the limits given 
by the subadditive ergodic theorem as well as the precise way in which the subadditive quantities 
control the solutions of (jl.ip . The former must necessarily be handled by probabilistic methods 
and the latter by pde methods. 

These are not mere technical difficulties. It turns out that, in the stochastic setting, the magni- 
tude of the fluctuations of the solutions u £ about their limit u must be separated into two separate 
regimes, which we refer to as ballistic and sub-ballistic, respectively (this terminology is borrowed 
from the probability literature, see, for example, Sznitman [34]). Intuitively, in the ballistic regime, 
the solutions are able to "feel" the random environment sufficiently quickly as e — > 0. Then the 
mixing of the medium dominates, which results in an algebraic convergence rate. In the sub-ballistic 
regime, the dependence of the solutions on H is highly localized in the vicinity of points y in which 
H(p,0,uj) is close to its maximum. It is therefore the law of H(p,0,-) near its maximum that 
principally governs the rate at which homogenization occurs, and this rate may be arbitrarily slow 
without a further assumption on the law. 

To give a more detailed overview of the approach, we start from the approximate cell problem 

Sv s + H(p + Dv s ,y,io) = in R d , (1.4) 

which, for each fixed p € M rf and 5 > 0, admits a unique bounded, uniformly continuous solution 
v 5 = v s (y , oo ; p) . The introduction of fll-4|) in the context of the homogenization of Hamilton-Jacobi 
equations goes back to the original proof of periodic homogenization due to Lions, Papanicolaou 
and Varadhan [21]. It is well-known by now (see, e.g., [3]) that the homogenization of (II. lj) to (II. 2p 
is equivalent to (and the effective Hamiltonian H can be identified by) the limit 

lim -Sv 5 (0,u ;p) = H(p) P-a.s. (1.5) 

Moreover, this equivalence is easy to quantify in the sense that an error estimate or convergence 
rate for (jl.5p can be transformed into one for homogenization. We are therefore left with the task 
of quantifying the limit in (jl.5p . 

The intuitive reason for the difficulty of arguing directly for the limit (|1.5p in the random case 
is the complicated dependence of the v s> s on H, which is both singular (information propagates 
only along characteristics and does not spread out) and global (information may travel far away 
in space, which is compounded by the lack of compactness). This problem is overcome in the 
stochastic setting by (i) imposing some kind of convexity assumption on H and (ii) using the 
subadditive structure of the metric problem (or its time-dependent analogue) to obtain an almost 
sure limit via the subadditive ergodic theorem. A comparison argument (introduced in [3], [4] ) then 
yields that, in the ballistic regime (p's satisfying H(p) > minif), the metric problem controls the 
limiting behavior of the 5v s (0,uj ;p) as 5 — > 0. In the sub-ballistic regime, i.e., for p's belonging to 
the "flat spot" {H(-) = mini/}, the limiting behavior of the 8v s is driven primarily by the law of 
H(p,0, •) near its maximum value, as mentioned above, and it turns out to be the thickness of the 
tail of this distribution which governs the rate of homogenization. 

We continue by introducing the metric problem: for each fixed /j, large than a certain constant, 
which turns out te be mini/, x € M. d and ui € there exists a unique nonnegative continuous 



-1 



S. N. ARMSTRONG, P. CARDALIAGUET, AND P. E. SOUGANIDIS 



solution = m At (-,x,a;) : M d — > R of 

H(Dm^, y, u) = fj, in R d \ {1} and m^(x,x,io) = 0. (1.6) 

In terms of control theory, the quantity m^y, x, u) corresponds to a "cost" of transporting a particle 
from x to y in the medium uj. It thus has the properties of a metric and is analogous to the time 
constant in first-passage percolation (see Remark |A.3|) . 

The m At (-, x, w)'s are the maximal subsolutions of H(Dw,y,u)) < fi in M. d subject to w(x) = 0, 
and this implies a subadditivity property. An easy application of the subadditive ergodic theorem 
(see [I]) then yields the existence of € C(M. d ) such that 

lim t _1 m M (iy,0,w) = mJy) P-a.s. (1.7) 

The deterministic rn^ can be identified as the unique nonnegative solution of 

H(Dm^) = ii in R d \ {0} and m M (0) = 0, 

and this provides another way of identifying H, one level set at a time. (See Appendix [A] for more.) 

The existence and some basic properties of the m At (-, x, cj)'s have been known for some time 
(see, for example, Lions |20j). More recently a simple comparison argument was introduced in [1] 
which demonstrated that the m At (-,x,w)'s control the 5v s (-, ui ;p)'s from below for every p G W d , 
and from above for p's in the ballistic regime. It follows from this analysis that the limit (jl.7|) 
implies the homogenization of (jl.ip . As we show here, this argument is constructive in the sense 
that a quantitative rate for the convergence of (I1.7P implies a rate for fjl .5f) . 

The main advantage of the metric problem is that it is localized. Indeed, while changes in the 
medium may influence the value of 5v s at far away points, the quantity m^(y,x,ui) depends only 
on the values of H(p,z,u) for z's satisfying m^(z, x,u) < m fM (y,x,uj), which is a bounded set with 
diameter proportional to \y — x\. This localization (see Lemma I A. 51 and (|A.28p below) permits us 
to use the independence of the medium by way of the martingale method of bounded differences 
and an application of Azuma's concentration inequality in a similar manner as in first-passage 
percolation [15] . 

In addition to quantifying the limit in fjl -Tp and hence in (|1.5p for p's in the ballistic regime 
(and obtaining an almost sure, algebraic rate of convergence), we identify a necessary and sufficient 
condition (see (|2.1ip below) for an algebraic rate of convergence to hold for p's in the sub-ballistic 
regime. This follows from a direct analysis of the v s, s using explicit comparison arguments. Merging 
the results for the ballistic and sub-ballistic regimes then yields, under assumption (|2.1ip . an 
algebraic rate for (jl.5p for all p's. 

We remark that we do not expect our arguments to yield sharp error estimates or convergence 
rates for homogenization. Indeed, as we explain in Remark 14. 2\ this is related to outstanding 
conjectures on the fluctuations of the time constant in first passage percolation. 

As mentioned above, the periodic homogenization of Hamilton-Jacobi equations was proved 
in [21]. This was simplified in [12] and subsequently extended to almost periodic media by Ishii [17J. 
The stochastic homogenization of convex first-order Hamilton-Jacobi equations was first proved 
in [321 ETJ and, for viscous convex Hamilton-Jacobi equations, by Lions and Souganidis [23J and 
Kosygina, Rezakhanlou, and Varadhan [19 1 . Lions and Souganidis [22] obtained results on the 
existence and nonexistence of correctors in the random setting and introduced in [21] a more direct 
proof of homogenization in probability. Later, more direct proofs of almost sure homogenization, 
based on the metric problem, were given in [31 0]. 
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The metric problem has been used by Davini and Siconolfi [10\ [TT] to study some connections 
between the stochastic homogenization of Hamilton-Jacobi equations and weak KAM theory and, 
for periodic H 's with special structure, by Oberman, Takei and Vladimirsky [29J and Luo, Yu and 
Zhao [25] in order to implement efficient numerical schemes for computing H. 

Outline of the paper. In the next section we give the precise assumptions and the statement 
of the main results. In Section [3] we review some preliminary results needed in our arguments. 
Controlling the fluctuations of the metric problem is the topic of Section [U These estimates are 
combined with comparison arguments in Section [5] to obtain corresponding bounds for the approx- 
imate cell problem in the ballistic regime. The sub-ballistic regime is studied in the second part of 
Section [5l where we produce error estimates under an auxiliary hypothesis on the law of H as well 
as examples demonstrating that, without such a hypothesis, the rate may be arbitrarily slow. We 
complete the proof of the error estimates in Section [6] and give convergence rate for the homogeniza- 
tion of the time-dependent problem (jl.ip . Finally, in Section [7] we discuss the convergence rates of 
the homogenization of (|1 .4[) and (jl.ip in almost periodic media. In the appendices we summarize 
the fundamentals of the metric and approximate cell problems. 

Notation and conventions. The symbols C and c denote positive constants which may vary 
from line to line and, unless otherwise indicated, depend only on the assumptions for H and other 
appropriate parameters (often an upper bound for \p\ or fi). For s,t E R, we write sAt:= min{s,i} 
and sVt := max{s, t}. We denote the d-dimensional Euclidean space by R d , Q d is the set of elements 
of R d with rational coordinates, N is the set of natural numbers and N* := N\ {0}. For each y £ R d , 
\y\ denotes the Euclidean length of y. If E C M. d , then \E\ is the Lebesgue measure of E, int E the 
interior of E, E the closure of E and conv E the closure of the convex hull of E. For r > 0, we set 
B(y,r) := {x € R d : \x — y\ < r} and B r := B(0,r). The distance between two subsets U, V C R d 
is dist(C/, V) = inf{|x — y\ : x € U, y € V}. If / : E — > R then we denote osce / := sup E f — infe /. 
If if is a finite set, then \K\ is the number of elements of K. The set of Lipschitz functions on 
a set U C R d is written Lip(C7) = C°' 1 (f7) and we set £ := Lip(M d ). The set of bounded and 
uniformly continuous real- valued functions on a metric space Y is denoted BUC(Y), and USC(Y) 
and LSC(y) are respectively the sets of real- valued upper and lower semicontinuous functions on Y . 
The Borel cr-field on R d is 23. If Si and S2 are cr-fields on sets X\ and X2, respectively, then Si<X'S2 
denotes the <r-field on X\ x X2 generated by 9i X 32- For a probability space (0,3", P), we say 
that an event A £ 3" is of full probability if PLA] = 1. We denote the indicator random variable of 
A G 3" by l A . If E G S <g> J and x € R d , then l {xeE} is the indicator random variable of the event 
{uj € £1 : (x,co) € E}. If X is a random variable and S C 3" is a cr-field, then E [X\S] denotes the 
conditional expectation of X with respect to 9- 

2. The assumptions and the statement of the main results 
We introduce our hypotheses and state the main results of the paper. 

2.1. The hypotheses. Let (f2,3~, P) be a probability space endowed with a group {T y ) yGK d of 3~- 
measurable, measure-preserving transformations r y : £1 — > £1. That is, we assume that, for every 
x, y G R d and A E J, 

f [ t j/(^)] = F I^] and (2.1) 
The Hamiltonian H : R d x R d x — > R is assumed to be measurable with respect to a S <g> "B ® 3". 
We write H = H(p, y, uj) and require that H be stationary in its dependence on (y, uj) with respect 
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to the translation group (r y ) y ^d, that is, we assume that, for every p,y, z G R d and wGfl, 

H(p,y,T z Lu) = H(p,y + z,uj). (2.2) 

In order to state the finite range of dependence assumption, which means roughly that H is "i.i.d." 
in its spatial dependence, we let 5(E) denote, for each E G "B <8> 3", the cr-field on Q generated by 
elements of 3~ of the form 

{oj G Q : H(p,x,oj) <a}n{uefl : (x,uj) G E 1 '} , 

where E' G 2 ® 3" such that E' C\ E, p,x £ R d and a G R. Intuitively, 9(E) consists of the 
information one learns by observing the Hamiltonian H(p, x, oj) for every p G R rf and x G {y G R d : 
(y, w) G -E}. If V G B, then we put 9(V) := 9(V" x fl) and we may henceforth assume, by restriction 
if necessary, that 3" = 9(R rf )- The finite range dependence hypothesis is then the requirement that 
there exists a positive D > such that, for all V, W G 23, 

if dist(V, W) > D then 9(F) and 9(W) are independent. (2.3) 
Of course, this implies that the group (T y ) y£Md is ergodic, but is much stronger. 

We continue with other structural hypotheses on H. We assume, for each R > 0, that the family 

{H(-, -,oj) :w G ft} is precompact in C(B R x R d ) (2.4) 

and 

{H(-,x,oj) : oj G 0, x G R d } is bounded in C ' 1 (Br). (2.5) 
We also require that -ff is coercive in p, uniformly in (y,oj), i.e., 

lim essinf H(p, y, oj) = +oo. (2-6) 

|p|->oo (i/,a;)e]R d xn 

We assume that i? is slightly more than level-set convex in p. Precisely, we assume that there 
exists A:lxK->l, which is nondecreasing in each variable, such that, for all (i,i/El, 

A(fj,,v) < fi\/ v and A(fj,,u) < fiV u if v ^ fj,, (2-7) 

and that H satisfies, for all p,q,y G R rf and wGfl, 

# + <A(H(p,y,u),H(q,y,u)). (2.8) 

Of course, is convex if and only if (12. 8j) holds with A(/i, f) = + z/). 

We also make the following assumptions regarding the shape of the level sets of H: for every 
p, y G R d and w G O, 

H(p,y,uj) > H(0,y,u) and esssupF(0, 0,w) = 0. (2.9) 

From the point of view of optimal control theory, the fact that there is a common po for all oj at which 
H(-, Q,oj) attains its minimum provides some "controllability", i.e., upper and lower bounds on the 
length of optimal paths. We loose no generality by assuming po = and esssup wer2 H(0, 0, oj) = 0. 
From our point of view, (|2.9p controls the growth of the m^'s (see (|A.10|) in Appendix |A|) . 

With the exception of Section [71 the hypotheses (|2.ip - (|2.9p described above are in force through- 
out the paper. For ease of reference, we write 

(I21|> . (gjj) , (|2~3l) . (jZaD , fl2J>D , ([23D , dZZD , flZSD and (J23]) hold. (2.10) 
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Some of our results are proved under an extra assumption on the distribution of H (0, 0, •) near 
its maximum. Precisely, this extra hypothesis is that there exist 9 > and c > such that, for 
every < A < c, 

P [H(0, 0, ■) > -A] > c\ e . (2.11) 

In light of (|2.9p . we see that, roughly speaking, (|2,lip is a requirement that the event that H(0, 0, •) 
is near its maximum is not too unlikely. For example, if H(0, 0, •) attains its maximum on a set of 
positive probability, then of course (|2.1ip holds for 6 = 0. 



2.2. The main results. Our first main result consists of error estimates for the limit (|1.7p . which 
measure the likelihood that the quantity {m^y, 0, w) — m^(y)\ is large relative to \y\. The definition 
and basic properties of the metric problem (jl.6p and its solutions and are reviewed in 
Appendix [Aj The proof of Theorem Q] is given in Section 

Theorem 1 (Error estimates for the metric problem). Assume (|2.10p . For each K > 0, there 
exists constants C\,C2 > 0, depending on H such that, for every 0<fi<K,X>0 and \y\ > D, 



"V(2/,0, •) -m^y) < -A 



< exp 



/xA 2 

cm 



and, if 



then 




A > C, 



log 1 + 



< exp 



/uA 2 



(2.12) 



(2.13) 



(2.14) 



The constant C\ is given in Proposition \4-l\ and C2 



4Ci\y 

2C14, where C14 is given in Lemma\4.10 



The error estimates for the metric problem and a careful quantification of the comparison ar- 
guments introduced in [2], together with an analysis of the convergence on the flat spot under 
the additional assumption ()2.1ip . yield the following error estimates for the limit (|1 . 5j) . The basic 
properties of the solutions v s of the approximate cell problem (jl.4p are outlined in Appendix [Bl 

Theorem 2 (Error estimates for the approximate cell problem). Assume (|2.10p . For each K > 0, 
there exists C3 > such that, for every \p\ < K, the following hold: 



C,5 



(i) For every < 8 < A < 1, 
-5v s (0, -;p)> H{p) + A] < C 3 5' 3d exp 

(ii) IfH{p) > and A, 5 > satisfy < 5 < 1 and 

A > C 3 (H( P )-§65 +H(p)- 1 5 1 ^ (1 + I logo) + \logH(p)\)\ 



then 



Sv d (0, • ; p) < H{p) - A < C 3 5~ 3d exp 



H(p)X 2 



(2.15) 

(2.16) 
(2.17) 
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(iii) Assume (|2,lip . There exists c\ > such that, if 5 > and < A < c\ satisfy 

A > C 3 ^|lo g( 5|i, (2.18) 

-5v s (0, -;p)< H(p) - A] < C 3 <5- M exp f-i- (j A . (2.19) 

By a covering argument and an application of the Borel-Cantelli lemma, the error estimates 
contained in Theorem [2] yield P-almost sure, local uniform rates of convergence for the limit (jl.5p . 

Theorem 3 (A convergence rate for the approximate cell problem). Assume (|2.10p . For all K > 0, 

there exists an event f2i G 3~ of full probability and a constant C\ > sitc/i i/iai ; /or every \p\ < K 
and U) & Ox, £/ie following hold: 



(i) For even/ i? > 0, 



hmsup sup vy ' — < 1. (2.20) 

<5^0 y&B R/s C4<5 3 | log (5| : 



3 



(ii) If H(p) > 0, then, for every R > 0, 



hminf inf " " ^ > -1. (2.21) 

*-+o j/eB B/5 C^Cp)- 1 ^ | log <J| 5 



(iii) // (I2TTT1) /joWs and 



a := - A - — - and /3 := - , (2.22) 
6 d+9 H 4 ' v ; 



i/ien, /or every i? > 0, 



hminf inf ,^ "f > -1. (2.23) 

The previous two results are proved in Section [5l where we also give a converse to Theorem [3|hi), 
which states that the extra assumption (|2.1ip is actually necessary for an algebraic rate of conver- 
gence to hold at p = 0. Indeed, keeping in mind that our assumptions imply that H(0) = 0, we 
prove in Proposition 15.71 that, roughly, if (|2.1ip is false, then, for every exponent rj > 0, 

1 . m . nf -^(0,^;0) = _ oo 

Furthermore, for any modulus function p, we construct examples of i/'s satisfying (|2.10|) for which 

hminf -^(°^;°) <_l F _ a . s . (2.25) 

It is therefore necessary to impose, in addition to (|2.10p . some assumption on the distribution of 
H (0, 0, ■) near its maximum to obtain a rate for the limit (|1.5p at p = 0. 

We next present our main quantitative results for the homogenization of (II. ip . Here u £ and 
u denote, respectively, the unique solutions of (11, ip and (jl.2j) subject to the initial condition 
u £ (-,0) = u(-,0) = no € £, which are bounded and Lipschitz continuous on M d x [0, T] for each 
T > 0. We begin with exponential estimates for the probability that \u e {x,t) — u(x,t)\ is large. 
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Theorem 4 (Error estimates for homogenization) . Assume (|2.10p . For each K > 0, there exists 
a constant C5 > such that, for every uq G C 0,1 (1R c( ) satisfying ||wo||c°» 1 (R d ) — K an d T > 1, the 
following hold: 



(i) For every < e < 1 and A > C§£3 , 

t\\ < —\T < n.T Gd \ 9d F- 6d P Y n 1 — 

C 5 e 



inf inf (u e (x, i, •) - u(aj, t)) < -AT 

xeB T 0<t<T 



< C 5 T M A 9fl! e~ M exp ( ) . (2.26) 



(ii) // (|2.1ip holds, then, for every < e < 1 and 

A > C 5 e^|loge|^, (2.27) 

sup sup {u £ (x, t, •) — u(x, t)) > XT 

x&B T 0<t<T 

< C7 5 T M A M e~ M exp ( -i- ( — A ) ) . (2.28) 



C 5 V e e d 

Our final main result is an almost sure, locally uniform, algebraic rate of convergence for the 
homogenization of (jl.ip . 

Theorem 5 (Convergence rate for homogenization). Assume (|2.10p and let K > 0. There exists an 
event Q2 G 3~ 0/ /nZZ probability and a constant Cq > suc/i i/iat, /or ei>en/ w € O2 a?id i*o £ C 0,1 (IR rf ) 
with \\uo\\c°' 1 (R d ) — K> the following hold: 

(i) For every T > 1, 

liminfinf inf ~ > -C 6 T. (2.29) 

e^O x€B T 0<t<T £5 I log e| 5 

(ii) If (|2,lip holds, a and (3 are as in (|2.22p and 

aid - /3 3 o! + 6» 

a:= = -A and b : = — - — = — V — — , (2.30) 

l + 2a 8 3d + 9 1 + 2a 16 4(3o! + 9) ' v ; 

then, for every T > 1, 

u e (x,t,uj) -u(x,t) 
limsup sup sup = < CeT. (2.31) 

e^o x&B T o<t<T e a |loge| fe 

Remark 2.1. We discuss later the sharpness of the exponent a and b. Let us point out for the 
moment that, in the special case that H is positively homogeneous of order one in p, i.e., for every 
t > 0, p,y G R d and u G Q, 

H{tp,y,u) = tH{p,y,uj), (2.32) 
the condition (|2.1ip is clearly satisfied for 9 = 0, and thus Theorem [S] gives a rate of Ofea | loge|i6) 
for homogenization. Moreover, this rate can be improved since f|2.32|) implies that H is also pos- 
itively homogeneous of order one, or equivalently that fi >— > m^y) is positively homogeneous of 
order one, that is, for all \i > 0, x, y G M. d and ui G ft, 

m^(y,x,uj) = fim 1 (y,x,uj). 

Thus the fluctuations of m^ — m^ are proportional to /i, and this prevents ()2.14p from degenerating 
as fi — >• 0. Indeed, we find that (I2.14p holds for every A > satisfying 

A >C 2t i\y I § (log(l + |y|)) * 



10 



S. N. ARMSTRONG, P. CARDALIAGUET, AND P. E. SOUGANIDIS 



instead of the more restrictive (|2. 13j) . One may propogate this improvement through the rest of the 
arguments in the paper to find, for example, that f|2.23|) holds for a = | and /3 = |, and (|2.3ip for 
a = ^ and b = ^j. The same observation holds for Hamiltonians which are positively homogeneous 
of any positive order, and we expect that other such improvements are possible for i?'s with special 
structure. 

3. Preliminaries 

We continue with some elementary remarks regarding our assumptions. First, it is easy to check 
that (|2,9p implies that, for each u) belonging to a set of full probability, 

sup H(Q,y,u) = 0. (3.1) 

Indeed, that the left side of (|3.ip is almost surely constant is a consequence of the ergodic hypothesis 
and the translation invariance of this quantity. That this constant must be zero follows from 
the continuity of H and the fact that (|2.9[) implies that P [sup yg Qd H(0, y, •) < 0] =1. Similar 
observations can be made for the sups and infs of other stationary random variables. 

Throughout the paper, all differential inequalities are taken to hold in the viscosity sense. Readers 
not familiar with the fundamentals of the theory of viscosity solutions should consult standard 
references such as E] . 

3.1. A discretization scheme. In the proof of Proposition 14.11 below, it is useful to employ a 
discretization scheme which for convenience we describe here. 

As we discuss in more detail in Appendix El the reachable set is defined for each fj,, t > and 
x € R d by 

&v,t( x ) : = {(V^) 6 M d x Q : m^(y,x,u}) < i} , 

and we also set 

^ t (x):={yeR d : (y,u)eK^(x)}, 31^(0) and B£ t := B£ t (0) (3.2) 

We also define, for each fi > 0, the filtration {9>,t} t>0 by GF^o := {^,0} and, for every t > 0, 

:= 9 (3W) • (3.3) 

That 1^^ is increasing in t follows from the fact that 9(E) is monotone with respect to E, according 
to its definition. 

We next describe the discretization scheme, which allows us to essentially condition on the 
identity of the reachable set Jl^t and apply the independence assumption (|2.3p in the form of 
Lemma |3.1| below. 

Define, for every r > 0, 

X r :={Ae-B : A = AQB r } 

and note that JC r is a compact metric space under the Hausdorff distance (c.f. Munkres |28j). which 
is defined by 

distn(E, F) := inf sup \x-y\V inf sup \x - y\ = inf {e > : E C F + B e and F C E + B £ } . 

x^EyeF y^ F x&E 

Note that, for every E,F G S, 

disttf (E, F) > dist (e, R d \ f) . (3.4) 
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Fix a small 5 > 0. Then, by the compactness of % r , there exists £ = 1(5, d,r) G N and a disjoint 
partition T\, . . . ,Tg C % r of 3C r into Borel subsets with diam#(Tj) < 5. Let Ki E 23 be the closure 
of the union of the elements of Tj. Then for each A G 3C r , there exists a unique 1 < i < £ such that 
A E Tj, which in turn implies that ^4 C Also define 



if* := Ki + 5 D 

so that dist(-ftTj, M. d \ Ki) = D . We have arranged things so that, for every 1 < i < 



S(Ki) and S R"\Jfj are independent 



(3.5) 



(3.6) 



and, for each A E % r and 1 < i < £, 

tiAeTi, then A C Ki C Ki C A + £d +< 5. 

We remark that the partition {I\} as well as the -fQ's depend on r and 5, but for convenience 
we do not explicitly display this dependence. 

The following lemma captures the intuitively obvious assertion that the behavior of the medium 
inside the set "Rati conditioned on the event that 31^ t E Tj (and hence 9l>u,t C Ki) is independent of 
the behavior of m fJi (y, Ki,cj), where we recall that the latter is defined in (|A.30|) and is independent 
of 9(Ki) by ()A.34p . Roughly, it contains a processed form of independence assumption which is 
particularly well-suited to our needs in the proof of Proposition 14.11 



Lemma 3.1. For each l<i<£,t>0 and A E 3>,t, 

1 An{X M , t 6r i } is S (Ki) -measurable. 

Moreover, 



E 



mf,(y,Ki,-)l { ^ ter . } 



E 



m^y^Ki,-) l {3i|titer .}. 



(3.7) 
(3.8) 



Proof. It suffices to show (|3,7p for A of the form 

A = {uj E : i?(p, x, w) < a} n {w G : (x, w) G -E} , (3.9) 

where E G 23 <8> 3" is such that .E7 C 3?^, p, x G R rf and a G R, since such events A generate 3>,t- 
Recall the definition of m^f in ()A.24p for nonempty and closed K C R rf and the fact from (|A.27p 
that, assuming G K, we have 

(y, 0,u)>t on dif implies that = |y G R d : m* (y, 0, w) < t| . (3.10) 

Due to (|32J| and (|3~1Q|) . 

in{we!l : ttj^ G Ti} 

= {a; G O : H(p,x,u>) < a , (x,w) G -E, {m^(-,w) < t} G Tj and m^(x,a;) < i} . 

If x then it is impossible that both (x, u) G E C and {m5' (■, w) < t} G Tj by ()3.10p . Thus 
this set is empty (and hence belongs to S(-fQ)) in the case that x €" Ki or {m^(-,w) < f} IV It 
is also S(-fQ)" measura t ) le if x E if j and {m^ l (-,uj) < t} G Tj because in this case £" C 3?^ )t C Ki xU 
and the map ui 1— >■ m Ki (x,uj) is S measurable. This confirms (|3.7p . 

According to f|3.5|) . (|A.34|) . and (j3.7|) . for each A G S^t, the event An {31^^ G Tj} is independent 
of the random variable m^iy, Ki, ■) . Hence for each A G 3>,t, 



E 



m^(y,Ki,-)l An{%ii teTi } =E m^(y,Ki,-) F [A n {tt„, t 6 . 
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The claim (13.81) now follows. □ 



3.2. Some preliminary results. To control the oscillations of solutions of the metric problem 
around their means, we use the "martingale method of bounded differences" based on the Azuma 
inequality pj]. See McDiarmind |27j and Alon and Spencer [2] for an overview of this probabilistic 
method, as well as a proof of Azuma's inequality, which is stated as follows. 

Proposition 3.2 (Azuma's inequality). Let {Xk}keN be a discrete martingale with Xq = 0. Assume 
that there exists a constant A > such that, for each k £ N ; 

ess sup |-Xfc+i — Xk\ < A. 
n 

Then, for each A > and N > 1, 

X 2 



¥[\X N \ > A] < exp 



2A 2 N 



We next state Hammersley's generalization of Fekete's lemma on subadditive functions. For a 
proof, see [15^ Theorem 2]. 

Lemma 3.3 (Hammersley-Fekete lemma). Suppose that £ > and f : [£, oo) — > R satisfies, for 
every s,t>£, 

f(s + t)>f(s) + f(t)-A(s + t), 
where A : [£, oo) — > R is nondecreasing such that 

"°° A(s) 



s 2 



■ ds < oo. 



Then r := lim^oo f(s)/s G (-co, oo] exists and, for every t > 

T >m + M_ 4 r^u. 



t ' t J 2t s 2 

Several of our arguments rely on the comparison principle for viscosity solutions of first-order 
equations, typically in the following form (see [9] or [S] for a proof). 

Proposition 3.4. Let G € C(R d x R d ), U be a bounded open subset ofR d , and u, -v € USC(ZT) 
and f £ C(U) satisfy 

G{Du, y) < f(x) < G(Dv, y) in U. 

Then 

supfii — v) = maxfu — v). 
u 8U 

We conclude with the following helpful lemma, which is due to the level set convexity of H and 
is useful for checking whether u € L is a subsolution of the equation H(Du, y, oj) < fj, for fx € R. A 
simple proof can be found in [1]. 

Lemma 3.5. Let \i € R, u) G Vl and U C R d be open. Then u G USC(£7) is a viscosity solution of 

H(Du,y,iu)<n inU (3.11) 
if and only if u is locally Lipschitz in U and satisfies ([3. lip almost everywhere in U. 

A commonly used fact in the theory of viscosity solutions is that the supremum (infimum) of 
a family of subsolutions (supersolutions) is a subsolution (supersolution), see |9j. Observe that, 
in light of Lemma 13 . 5 (. the infimum of a family of subsolutions of (|3.1ip is a subsolution and, in 
particular, the infimum of a family of solutions of (13. 5^ is a solution. 
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4. Error estimates for the metric problem 



Here we study the fluctuations of about their homogenized limit m„. The ultimate goal is to 
obtain the exponential error estimates contained in the statement of Theorem [TJ We assume that 
H satisfies (|2.10p . but we do not assume (|2.1ip . We also fix K > 1 and 

< n < K. (4.1) 

The positive constants in our estimates denoted by C, C\, C*2, . . . and c, c%, C2, . . . depend only on K, 
the underlying dimension d and the assumptions for H. Several of our estimates depend on a lower 
bound for fi, and since we must keep track of this dependency, we explicitly display dependence 
on [i. We allow the constants C and c to vary from line to line. 

There are essentially two steps in the proof of Theorem [TJ The first is to obtain exponential 
error estimates controlling the oscillations of m^(y,0, •) about its mean 

M^y) :=E[m M (y, (),•)]. 

The second step, which is more involved, is to estimate the difference between the deterministic 
quantities M^(y) and m^(y). 



4.1. Controlling the oscillations around the means. We begin by estimating, for large \y\, 
the probability that m^y, 0, •) is far from its mean. We use an argument inspired by the pioneering 
work of Kesten [18] in the theory of first-passage percolation, who introduced a martingale method 
based on Azuma's concentration inequality. We also benefit with some very elegant simplifications 
of the argument due recently to Zhang [37] , 

Notice that, unlike in percolation theory (or its continuum analogue), our Hamiltonian is not 
assumed to be positively homogeneous. In this generality, it is necessary to keep track of the 
dependence of the estimates on a lower bound for (jl. We recall that, in view of Proposition IA.2f iv) 
and ()4.ip . there exist constants 

< c/i < 1^ < Lf, < C (4.2) 

such that, for every x, y £ 

ly\y ~ x\ < m^(y,x,(jj) < L^\y - x\. (4.3) 

In the control theory interpretation (see Remark lA.3p . this important estimate, which we use many 
times below, provides upper and lower bounds on the lengths of optimal paths connecting two 
points x, y £ M d . 

Proposition 4.1. For each A > and \y\ > D, 

K(i/,0,-)-M M (i/)|>A] <exp(-^|), (4.4) 

where C\ > is given explicitly by C\ := 8^Z~ 1 -L^(2Z /J + L^D < 24fj,l^ 1 L^ l D < CL^D. 

Proof. Fix y £ M. d with \y\ > D and define 

T:=L M |y| and r := T/l^ = (4.5) 
so that (see (|A.18[) ). for every u £ J7, 

$%,t £ B r . (4.6) 
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We employ the discretization scheme, for fixed 5 > 0, as described in Subsection 13.11 with the 
notation introduced there, and we define a continuous-time martingale {Xt}t>o adapted to {9"^,t}t>o 
by setting, for each t > 0, 

X t :=E[m^y,0,-)\^,t] - M^y). (4.7) 
Due to 3^,0 = {0, O}, T = L^\y\ and (EQ9|) . 

X = and X t (u) = m M (y, 0, u) - M^y) for every t > T. (4.8) 

Our goal is to apply Azuma's inequality in order to estimate the oscillations of Xt- We must first 
obtain an estimate of the form 

ess sup \X a (u) — X t {oj)\ < A + B\s — t\. 

We begin with the observation that, owing to (|A.29p . for every < t < s, 

E[m M (y,0,-)l{ tf 63i M }|3>,«] = m i*(yA-)1{v&li l >,t} ( 4 - 9 ) 

and, hence, 

X s - X t = E[m fl (y,0,-)t {y ^ t} \J llyS ] - E [m^y, 0, -)l {y ^ t} \ 3^*] ■ (4.10) 
Using (|A~20|) . we find that 

m M (y, 0,u)l [y ^ t} (u)) = min (m M (y, z,w) + m^z, 0, u)) l {yjSR 

zecU V,t 

= [m^y^^u)) + t) l{ y ^3i Mit }H =m M (2/,^,o;)+tl{ 2/ ^, t }H 
and, since {y ^ 3^,t} € 3"^,t, we may simplify (I4.10j) to write 

X s - X t = E[m IM (y,Ji fMj t,-)\^,s] - E[m M (y, •) | 3>, t ] . (4.11) 
Clearly for every < f < s, 

mix(y,^„s,^) < mn(y,0^ t ,u) < (s - t) +m IM (y,J^ a ,u). 
Combining the last two lines, we obtain 



\X S -X t \ < (s - t) + |E [m M (y, •) | 3^ s ] - E [m M (y, 3^, t , •) | 3>, t ] 



(4.12) 



We next use the discretization scheme to estimate E [m^y, ^,t, •) | 3>,t] by approximating the 
integral represented by the expectation as a sum of characteristic functions. With Ki, Ki and Tj 
as described there, observe that, by (I3.6p . 



m M (y, ^, •Jl^.teTi} < ^,t, -^{S^ter,} < (^C + *) + nty(l/, 0) 1 {K f .,ter i }. ( 4 -!3) 

Taking the conditional expectation of (I4.13P with respect to 3~ Mi t and applying (13.8H . we get 



E 



< (l m (D + 5) + E [m^y, «i, -)J J l{^, te r l} (4.14) 
Since {Tj} is a disjoint partition of JC r , we also have 



E 



m^y,^.^, ■) 



Yl E [ m /*(f» 3l M,*>-) :IL {* M ,ter i 



1 



(4.15) 
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Multiplying (|4.14p by l{x Ms er } an d summing over the indices i and j yields, in light of ()4.15j) . 



15 



< E 



m^y, K^t, •) | S^t] - E ^> •)] l{^, te r i } 1 {^, s er J } < + ( 4 - 16 ) 



In the same way, after interchanging s for t and j for i, we also obtain 



< E [m^y, %„, s , •) | ? M , S ] - [m M (y, •)] H^, t er i} ^, s er j} < L^(D + 5). (4.17) 

It follows that 

|E[m„(y, 3^, s , •) | 9>,J - E[m^(y, J^ t , •) | J^ t ] | 

< L^(D + | E K^' ^' •) " Kj> •)] I ^^Werjl^er,}- (4.18) 

If, for some i, j = 1, . . . , £, there exists ui belonging to the event that 3l^ t G Tj and 3?^ )S G Tj, then 
disttf (k u K 3 ) < dist H {K u Kj) < dist H (tt£ t , tt£ a ) +25 < + 25. 



Using pll . we conclude that, for every i,j = !,...,£, 



E 



L, 



m^iy.Ki,-) -E m^(y,Kj,-) t { ^ teTi} t { ^ seT]} <-^(s-t) + 2L ll d. 



Combining this with (|4.18p and sending 5 — > yields 

\E[m^y, X^s, •) | 3>, a ] - E[m M (y, tt M , t , •) | 9^] | < L M D + ^(s - t). 



Finally, from (I4.12p . we get 



\X t -XJ <L U D + 



+ 1 )(s-t). 



(4.19) 



(4.20) 



Define a discrete martingale sequence Xfc := Xhk with /i := I / '(I M + L M ) and observe that, 
according to (|4.20p . for all € N, 

— Xk\ < 2L fl D. 

An application of Azuma's inequality (Proposition 13.2) ) yields, for every A > and N € N, 

— 1 

8LlD 2 N J 



X N \ > A 



< 



exp 



(4.21) 



Let N be the smallest integer larger than T/h so that Xn = X? = 77^(^,0, •) — M^(y). It follows 
that, since \y\ > D and T = L^\y\, 



h 

From (USD, (1^2Tj) and (14T22D we deduce 
|m M (y,0,-) -M„(y)| > A 



\X N \ > A 



< exp 



-A 2 /, 



81^(2^ + L M )Z%| 



(4.22) 



□ 
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Remark 4.2. Integrating ()4.4p yields 

poo _ 1 - POO / „\ \ £1 

var (m^y, 0, •)) = / P K(y, 0, ■) - M M (y)| > A^ < / exp --^ UA = — |y|, 
Jo L J Jo V °i 12/1 / 

which mirrors the bound on the variance of the time constant obtained by Kesten [18J for first 
passage percolation. Obtaining an optimal estimate for the fluctuations of the latter is a well- 
known open problem. It is conjectured that the variance of the time constant should behave, in 
dimension d = 2, like Of|y|3j for large |y|, and it is believed that the oscillations should decrease in 
higher dimensions. Nevertheless, it is still open in every dimension d>2 whether, for some a < 1, 
this quantity is bounded by 0(|y| a ) as |y| — > oo. We expect that it will be similarly challenging 
to prove such a bound for our quantity vax(m /i (y 3 0, •)), and still more difficult to find the optimal 
exponent for the algebraic rate of homogenization of (jl.ip . 

In analogy with the best known variance bound in first-passage percolation, due to Benjamini, 
Kalai and Schramm [7], we expect that an estimate of the form 

var (m M (y, 0, •)) < - (r*Tl) (4.23) 

can be proved, in dimensions d > 2, by an application of Talagrand's concentration inequality |35j. 
In fact, as we were completing the writing of this paper, we received a new preprint by Matic and 
Nolen [26] who have obtained, in a slightly different setting, a bound like f|4.23j) for a certain class 
of Hamilton-Jacobi equations in special i.i.d. environments. 

4.2. The strategy for estimating M„(y) — m„(y). Having estimated the oscillations of m^(y, 0, •) 
about its mean M^(y) in Proposition 14.11 in order to prove Theorem [1] it remains to estimate the 
rate at which the means t^M^ity) converge, as t — > oo, to their limit m^(y). On one side our task 
is trivial. Indeed, owing to (jXij) and (|X8|) . we have 

M M (y + z) = E [m^y + z, 0, •)] < E [m M (y, 0, •)] + E K(y + z, y, •)] = M^y) + M^z). (4.24) 

It follows from Fekete's lemma (Lemma 13.31 in the special case that A = 0) and (|A.36P that, for 
every y € M d , 

M At (y) > Mt- l M^ty) = lim t^M^ty) = m M (y). (4.25) 

t>l t— >oo 

The estimate (I2.12P is then immediate: 

Proof of (|2TT2"j) . The inequality follows from (|4T23|) and □ 

In order to prove (|2.14p we are confronted with the more difficult task of finding good upper 
bounds for M^(y) — m^(y), which is the focus of the rest of this section. A typical argument for 
such an estimate (see the discussion in Hammersley [16]) and the strategy we use here involves 
approximating M„(y) by another quantity which is superadditive. Fekete's lemma (or, more pre- 
cisely, Lemma \3.3\i is then applied "from the other side" to obtain an estimate on the deviation of 
this approximate quantity from its asymptotic limit. The desired estimate in terms of the original 
quantity then follows, depending on the quality of the approximation. This strategy was used by 
Alexander [lj in the context of first-passage percolation to obtain estimates on the deviation of the 
expected passage time from the limiting time constant. 

In our context, it turns out to be more convenient to first obtain estimates for the difference 
between the quantities 



E[m„(H t ,0,-)] :=E 



min irin(z, 0, 



and rrin{Ht) '■= min m„(z) 
zeH t 
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where Ht is a given plane a distance t from the origin. We then argue that the value of E [m^(Ht, 0, •)] 
must be close to 

M^Ht) := min M^z), 
z&Ht 

which yields good estimates for the deviation of the latter quantity from m^(Ht). This is then 
transformed, using a simple geometric argument, into an estimate for M^{y) — m^(y) for large \y\. 

4.3. Introduction of the approximating quantity. Fix a unit direction e G M. d which for 
notational convenience we take to be e = e<i = (0, . . . , 0, 1). For each t > 0, define the plane 

H t :=te + {e} ± = {(xr,t) :i'£R i4 } (4.26) 

and its discrete analogue 

H t := j(n,t) : n G Z d_1 } . (4.27) 
We also denote, for t > 0, the halfspaces 

H+ = {(a/, x d ) G R d : x d > t J and JEr t ~ = j(x, s d ) G M d : x d < t} . (4.28) 
Define, for each <r, i > 0, the quantities 

G„, a {t) := V E [exp (-crm^y, (),•))] and ^(t) := -- log G M , CT (t). (4.29) 
yeH t 

We show that <? Mi0 -(i) is a good approximation of E [m M (i?t, 0, ■)] for suitable choices of cr > 0. We 
first need a technical lemma, also used many times below, which asserts that a substantial portion 
of the quantity G^it) is contributed by lattice points (n,t) G H t with \n\ < 0(t). This implies in 
particular that G^ a and g^ a are finite. 

Lemma 4.3. There exists C7 > such that, for every t>0,0<a<l and R > (L^/l^t, 

G^{t)<-^ Yl E[exp(-(7m M (y,0,-))]. (4.30) 
y eH t nB R 

Proof. According to (|4.3p . for every u G and y G M d , 

exp (-L^a\y\) < exp (-am^(y, 0, w)) < exp (— Z M <r|?/|) . (4.31) 

Thus 

^2 exp(-crm M (y,0,w)) < ^ exp (-l^a\y\) < ^ exp (-Z M cr|j/|) 

y&H t \B R y£Ht\B R y& d -^\B R 

<C r d - 2 exp(-l^ar) dr = Ca 1 ^ / exp (-i„r) dr, (4.32) 

J R JaR 

where C > depends only on d. Moreover, 

roo 

exp (— am^te, 0, w)) > exp (— L^at) = / exp(— l^r) dr. (4.33) 

Since R > (L^/l^t, it follows from (|4.32p and (|4.33p that, for every w £ H, 

exp (— am^(y, 0, w)) < Cer 1-d exp (— am^(te, 0, £j)) < Ccr 1_a! exp (— am fl (y, 0, cj)) . 

yeSt\B R yemnB R 
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Taking expectations yields 

G„At)= Y E[exp(-am M (y,0,-))]+ Y E [exp (-am M (y, 0, •))] 
yeH t DB R yeH t \B R 

<(l + Ca x - d ) Y E [exp (-<rm M (y, (),•))] ■ 
y eH t nB R 

We obtain flOO!) for C 7 := (1 + C), since <r < 1 implies 1 + Ca l ~ d < (1 + CV^. 



□ 



Next we show that g^it) gives a good upper bound for E [m^(H t , 0, •)] for large t and appropriate 
choices of a > 0. 

Lemma 4.4. There exists a constant Cg > smc/i £/ia£ ; for every t > D and < a < 1, 

E[m M (^,0,-)]-C 8 ^ + ilog(l + -^-^ <^(t)<E[m M (fl- t ,0,-)] + C 8 . (4.34) 

Proof. The upper bound in (|4.34p is easy. Using (jA.ll|) . we have 



GWt) > E 



sup exp(-crm M (y,0, •)) 



E 



> E 



exp (-cr (mJH t ,Q, •) + - 1 



exp -a inf m^(y,0, 



exp (-ffL„(d- l)^ E [exp (—amJH t , 0, •))] . 



After taking the logarithm of both sides of this inequality, an application of Jensen's inequality and 
a rearrangement yield the second inequality of (|4.34|) with Cs = (d — 1) 3 L„. 

To obtain the lower bound, we use both (|4.4[) and (|4.30|) . For every \y\ > D, we have 

E [exp (-am^y, 0, •)] = / cr exp(-as)P [m M (y, 0, •) < s] (is 

/■M M (j/) 

< cxexp(-o-M M (y)) + / a exp(-as)P [m M (y, 0, •) < s] 



A/ M (J/) 



1 + 



Applying (|4.4p . we obtain 



exp((rA)P[mu(i/,0,0 - M„(y) < -A] c?A a exp (-crMjy)) . 



/ f^fe) / uA 2 A \ 
E[exp(-(rm M (y,0,-)] < U + J exp (,A - ^ j dA I a exp (-aM^y)) 



We estimate the integrand above by completing the square, i.e., 

V 2 .. / L.l\ 2 



<tA 



Cih/I CiIj/I 



/x /■ x crC\\y\ V . 1 o „ , , , 1 



A 



2fi 



+ — ^CM < —a 2 CM 



and thus obtain 



E[exp(-crm M (y,0,-)] < ( 1 + M^y) exp ( — a 2 C x \y\ ) ) a exp {-aM^y)) . (4.35) 
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Summing ()4.35j) over y 6 H t PI Br, with i? := (L^/l^t, and applying (|4.30p . we get 



G^(t) < Cja 2 - d i 1 + M M(y)exp (^C'lll/lj) ex P (~° M M) 

y€H t nB R 

<C 7 a 2 ~ d (l + M M (y)ex P (i-a 2 C 1 |y|^exp(-aE[m / ,(^,0,-)]) 



y eH t nB R 

< C 7 a 2 - d R d - 1 exp (-<tE [m M (iJ t , 0, •)]) f 1 + exp " 

Since, in view of (|4.2p . /?, < Ct/fj,, we obtain 

G^(i) < CtV-V 1-d exp f-aE K(^,0, •)] + ^ 

Taking logarithms, dividing by — a and rearranging the expression yields the desired result. □ 

4.4. The (almost) superadditivity of <? Mj(J and estimates for E [m^(Ht, 0, u>)] — m^Ht). The 
next step is to prove that g^o- is essentially superadditive, which is summarized in the following 
lemma. Unlike the approach taken in [I], we do not use an abstract result like the van den Berg- 
Kesten inequality, which does not seem to easily apply in the continuous setting. We opt instead 
for a simpler "splitting technique" to apply the independence assumption more directly. A similar 
technique was employed by Sznitman [33]. 

The critical property of the m^'s needed here, which allows us to exploit the independence of the 
random medium, is the dynamic programming principle. It asserts that, if every path from x to y 
passes through a surface, then, for some z on the surface, the cost of moving from x to y is equal 
to the sum of the cost of moving from x to z and from z to y. Precisely, for every open f/Ci' ! 
with x £ U and every y € M. d \ U and well, 

mJy, x, u) = min (mJy, z, u) + mJz, x, u)) . (4.36) 

z&dU 

This is proved in Appendix [A] see Proposition IA.2f vi). 

Lemma 4.5. There is a constant Cg > such that, for every s,t> D and < a < 1, 

9,At + s)> 9^it) + 9»,*(s) ~ Q U + log (l + ^jf) ) . (4.37) 



Proof. Fix s, t > 0, y € H s+t , uj £ Q and observe that, in view of (I4.36|) 
m^{y,Q. 

Thus, using ([ATTH 



m^y,0,u)) = min (mJz,0,ui) + mJy, z,u>)) > mJH t ,0,u) +mJy,H t ,u)). 
zeH t 



m^(y,0,uj) > m^(H t ,0,u) + m^(y, H t+D ,uj) - L^D. 
Since € and y G , we have 

m^(H t ,0,uj) = m^(H^,Q,uj) and m M (y, H t+D ,u>) = m^y, H~ +D ,u), 
which is immediate from the expression (|A.30p and (|4.3p . Applying ()A.34p . we conclude that 

m^H^ ,0, ■) is 9(-fft _ )-measurable and m^(y,H^ +D ,-) is S(-fT t + ( _ D )-measurable. (4.38) 
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In light of (|2.3p . these random variables are independent and thus 

E [exp {-am^y, 0, •))] < exp (aL^D) E [exp (-am^Hf, 0, •))] E [exp (-am^y, H~ +D , •))] . 

Returning to the discrete setting, we next claim that, for R := L^(s + 1)/^, 

m M (i? t + ,0,u;) =m AI (-Hi,0,w) > min m^(z,0,u) - L^d - 1)5. (4.39) 

Indeed, it is clear from (|4.3[) that any z £ Ht attaining the minimum on the left side of (|4.39p must 
belong to Br, and ()4.39p then follows from (jA.lip . In a similar way, since \y\ < R, 

m M (|/,E ( " +fl ,w) > min m^y, z, u) — (d + (d — 1)5 

z£H t nB 2R v 

Combining these inequalities, we obtain 

E [exp (-0-771,^(^,0, •))] < exp(Ca) E [exp (—am^(z, 0, •))] E [exp (-am^y, z', •))] . 

2;,2;'G75 t nB2B 

Note that, if z 1 £ i^, then y — z' £ H s . So, in view of the definition of and (]A.4p . we have 

^2 E [exp (-am^y, z' , •))] = E [exp [—am^y — z', 0, ■))] = G MjCT (s). 
«'eH t z'&H t 

Therefore 

E [exp (-o-m M (y, 0, •))] < exp (Ca) G^(t)G^(s). 
Summing over all y € H t + S H Br and using Lemma 14.31 yields 

G^(s + t)< CR^a 1 ^ exp (Ca) G^(t)G^(s). 
Taking logarithms and dividing by a concludes the proof. □ 

We next use Lemma l3.3l to obtain a rate of convergence for the means t -1 E [m^(Ht, 0, ■)] to their 
limit m^Ht). 

Lemma 4.6. There exists a constant C\q > such that, for every t > D, 

E [m^Ht, 0, •)] <m^H t ) + C w ^log^l + ^V. (4.40) 



Proof. According to Lemma 14.51 t ne quantity g a ^ is almost superadditive. More precisely, for all 
s, t > 0, we have 

S<x,/*(s + t) > ^(s) + ft^i) - A a ^(s + t), (4.41) 

where, for Gg as in (|4.37p . 

A^i) :=^fl + logfl + — 



Since A CTi/x is increasing on [l,oo) and 

r°° a_..iTi 

■ dt < oo 



A g ,„(t) 

we may apply Lemma [3731 to deduce that g CT „ := limt_Kx> 9a,n(t)/t exists and, for every t > 1, 



1 J2t s 
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An easy integration by parts yields 



2/ 



s A at \ \ crfi 



In view of the second inequality in (|4.34j) , we also have 



< lim - (E [m^Ht, 0, •)] + CL^) = m^H x ) 
Combining the previous three estimates, we obtain 

< rn^Hx) + — f 1 + log ( 1 + — 
t at \ \ a/j, 

Multiplying by t, applying the first inequality in (|4.34j) and using the positive homogeneity of 
yields 

'at 1 1 . / t 



E[ mfl (H t ,0,-)] <m^H t ) + C — + - + -log 1 + — 

\ / u z a a \ 07/, 

and choosing a := fit~2 (log(l + i//i)) 5 completes the proof. □ 

4.5. Error estimates for Mu(y) — m^{y) and the proof of (|2.14p . It is the rate of convergence 
of t^ 1 M^ty) to m^(y) that we wish to estimate, not that of i _1 E [m^Ht, 0, ■)] to m^Hi). In 
order to reach our desired goal, we must compare the quantities M^(y) and E [m^(^,0, •)]. This 
is accomplished in two steps. The first is to show that the quantities E [m^^Ht, 0, ■)] and 

M^Ht) := minM M (y) 

are close, which then gives an estimate for the difference between M^Ht) and m^(Ht). The second 
step is to use elementary convex geometry to relate M^y) to the values of M^(H) for all the 
possible planes H passing through y. 

Lemma 4.7. There is a constant C\\ > such that, for each t > D, 

M^Ht) < E [m M (i?t, 0, •)] + Cn QlogTl + ^V. (4.42) 

Proof. Let R := (L^/l^t. For every uj € $7, there exists z & H t H Br such that m^(z, 0,uj) = 
m^(H t , 0, oj). Hence there exists £ € H t H Br such that 

171^,0, u) < m^(H u Q,uj) + L M (d -1)2. 
For every z £ H t we have E [m^z, 0, •)] = M^(z) > M^Ht) and thus, for every A > 0, 

{w€ft : M M (ff t ) -m M (# t ,0,w) > X + L^d-l)^ 

C |J <E SI : m fM (z, 0, w) < M^(z) — A} . (4.43) 

Applying (|4.4p . we find 

M^flt) - m^Ht, 0, •) > A + L^d - 1)^1 < CR d ~ l maxP [m^(z, 0, •) - M^z) < -A] 

J z£H t 

< CR d -' exp < C^H^ exp (-^) . (4.44) 
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We wish to use the expression 

roc 

M^H t ) - E [m^Ht, 0, •)] < / P [M^Ht) - m^H u Q, •) > A] 

Jo 

and then apply f[4.44j) to the right side of (|4.45p . but due to the factor t d ~ l on the right side 
of (|4.44p . this bound is not very helpful unless A is large relative to t. With this in mind we fix 
A > 0, to be selected below, define 



dX 



(4.45) 



Ai :=1 + 

and then estimate the right side of (|4,45p by 

/ P [M^Ht] - m^Ht, 0, •) > A] dX 
Jo 

< Ai + L M (d-l)5 + 



At 
/i 



log 1 + 



1-d.d-l 



< Ai + L M (d-l)5 + C^- a t 



M^Ht) - m^H u 0, •) > A + L^(d - 1) 
/xA 2 



dX 



cxp 



Ct 



dX. 



Observe that 



exp 



' Ct 



dx < ^- d t d - 1 



exp 



[iX\X 
Ct 



C^i-^-iJ-exp 
liXi 



dX 
Ct 



< c 



1 + 



A 

c 



By selecting A to be a large enough constant, the last expression on the right is at most C. 
Combining the last two sets of inequalities with (|4.45p . we obtain 



M^Ht) - E [m^Ht, 0, •)] < Ai + L^d - 1) 5 + C < Ai + C, 
which implies f[4.42j) . 

Lemmas 14.61 and 14.71 give an estimate on the difference of M^(Ht) and m^(Ht) 
Corollary 4.8. There exists C\i > such that, for every t > D, 



(4.46) 
□ 



M^H t )<m^H t ) + C 



12 



log 1 + 



(4.47) 



The relationship between M^{Ht) and M^{y) depends on the following geometric lemma. 
Lemma 4.9. There exists C13 > such that, for every N S N and a > 0, 

conv {y e R d : M^y) < a} C \y € M d : M^{Ny) < (N + C 13 //x)a} . (4.48) 

Proof. Let a > 0, G N* and select y € conv{y G M d : M^{y) < a}. According to Caratheodory's 
theorem (see for example [H]), there exist y\, . . . , yd+i and Ai, . . . , A^+i G [0, 1] such that 

d+l d+1 

y = ^2 ^jVjj 1 = XT ^ anc ^ M^(yj) < a for every 1 < j < d + 1. 
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For each 1 < j < d + 1, select aj £ Q such that o~jN € N and < Xj — o~j < 1/N, define 
z := X^jii a jVj an d observe that, by (|4.3p . 

i / d + 1 . . (d + l)a 

\v — z\ < max < 

- N i<j<N Wjl ~ Nip 

and, by (JO) and 



M M (iVy) < M M (7Vz) + L M iV|z - y| < M M (iVz) + ( v ' ? 7 g ) a < M^Nz) + — . ( 1.19) 



(d+l)L u \ ,„ T , Ca 



Notice that, since Mr,- € N and Nz = E^iC^jOVjj 

we may apply (|4.24j) to deduce that 

M^Nz) < ^(Na^M^yj) < Na^ °j < Na Yj X i = Na ~ ( 45 °) 
j=i j=i j=i 

Combining (|4.49j) and (|4.50p yields the lemma. □ 

The previous lemma and (|4.47|) yield a rate of convergence for M M (y) to m M (y). 
Lemma 4.10. There exists a constant C14 > such that, for every \y\ > D, 

M„(y) < m M (y) + C U ( + —) flog U + ^) V ■ (4.51) 



\ /i2 M / V V I* . 

Proof. The first step is to show that, for every z € M. d such that \z\ > D 



^econvjy Gl d : M„(y) < m^z) + C 12 l og (\ + Jfi^ 2 1 (4.52) 

where C12 > is as in (|4.47|) . Suppose on the contrary that (|4.52|) fails for some z € M. d with 
t := \z\ > D. By elementary convex separation, there exists a plane H with z £ H such that 

1 

t ( t \ \ 2 
M^(H) > m^z) + i4 where A := C u {-^ log ( 1 + - 

Since if is at most a distance of \z\ = t from the origin, we may assume with no loss of generality 
that H = H s for some s < t. We deduce that 

1 

M^H S ) > m^z) +A> m^Hs) + Ci 2 14 lo s ( 1 + - ' ' 

a contradiction to (fOT|) . Thus P~52]) holds. 

Now we fix y G ]R rf with |y| > D and apply (|4.52p to z := y/iV, with N € N* to be chosen below, 
to conclude that 



y/iV G conv M„{x) < m^y/N) + C 12 log ( 1 



2-1 
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and, after an application of (|4.48p . 



C 



M»{y) < I N + — ) I m^y/N) + C 



■ 12 



n fi- 



ll/ \ 



N\y\ 



< m^y) + Ci 3 L^ + C 12 I log ( 1 + ^- ) ) + (',,(',. 



■log ( 1 + 

w 



■log 1 + 



Nfi 



In the case that \y\~ < fJ, , then we let N be the smallest integer larger than \y\3 to get, after 
using that \y\ > D and \x < K, 

, 2 



M M (y) < m^y) + C 



\y\ 



log 1 + 



(4.53) 



If, on the other hand, < /i 3 < \y\ , then we let N be the smallest integer larger than \i 1 to get 

|y|a 



M M (y) < m M (y) + (log (1 + \y\)) 

f!2 



In either case we get (|4,5ip . 



(4.54) 
□ 



The previous lemma provides the desired estimate for the difference between t 1 M )1 (ty) and 
rn^{y) for large t > 0, and now the proof of Theorem [T] follows immediately. 

Proof of (|2.14p . The estimate is immediate from Proposition 14.11 and Lemma 14.101 □ 

4.6. Some further error estimates. We conclude this section with versions of (|2.12p and (|2.14j) 
which hold uniformly for y € Br. These estimates, which are needed in the next section, follow 
from Theorem Q] and a simple covering argument. 



Lemma 4.11. There exists C15 > such that, for every A > 4L„Z) and R > 3D, 



inf (m M (y,0, ■) - m M (y)) < -A 
yeB R 



< Ci 5 R d exp 



/1A 2 
AC X R 



and, if 



then 



A > 2Co 



R^ R,l 



3 



(J 



log 1 + 



sup (m M (y, 0, •) - m^(y)) > A 
yeB R 



R 
/' 



< C 15 R d exp 



/xA 2 
I6C1 



(4.55) 
(4.56) 
(4.57) 



Proof. We may select y 1 ,...,y N £ B R \B D with N < C(R/D) d < CR d such that B R is covered 
by the balls B(yj, 2D). Then by (lA.llj) we have, for any A > 0, 

N 



inf (m M (|/,0, •) - m M (y)) < ~ - 2L^D 
y£B R 1 



According to (j2TT2|) . for each 1 < j < N, 



A 



m M (y i ,0,-) -m^yj) < 



< 



exp 



/uA 2 



4Ci|l/j 



< exp 



/iA 2 
4Cii? 
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Since N < CR d , we obtain, for every A > AL^D, 

inf (m^d/,0,-) -m^y)) < -X 
y£B R 



<CR d e W < MA 



4Cii? 



The estimate (I4.57P is obtained in a very similar way from Theorem [T] and a covering argument. 
We omit the proof. □ 

5. Error estimates for the approximate cell problem 

Here we obtain estimates on the difference between —5v s (y,uj;p) and H(p), study the rate for 
the almost sure convergence 

lim sup 5v s (y, uj ;p) + H(p) = (5-1) 

<5^o y eB R/s 

and prove Theorems [2] and [3j The main difficulty comes from the fact that the rate for the 
approximate cell problem may be very different depending on whether p belongs to the flat spot 
{H = 0} or not. Recall that the flat spot is never empty since, e.g., H(0) = (see Appendix |A|) , 
Moreover, unless for example H is positively homogeneous with respect to p, the flat spot {H = 0} 
is in general equal to {0} (see, e.g., [3]). We use the metric problem to control the — Sv Sl s from 
above, and from below for p : s away from the flat spot {H = 0}. To obtain the upper bound on the 
flat spot we study directly the behavior of the Sv Sl s. 

The definition and basic properties of the v s (y, uj ;p)'s are given in Appendix [Bl We recall here 
two important estimates, namely 

— ess sup H(p, 0, •) < Sv s (y, uj ;p) < — ess inf H (p, 0, ■) (5-2) 
n n 

and 

v S (y,uj;p) - v s {z,u);p) <K p \y-z\, (5.3) 

where K p > is given explicitly in (|B,6h and depends only on H and an upper bound for \p\. 

Throughout this section the positive constants C, C±, C2, ■ ■ ■ and c,ci,C2, . . . depend only on an 
upper bound for \p\, the underlying dimension d and the assumptions for H. We allow the constants 
C and c to vary from line to line. 

5.1. The ballistic regime. We combine the exponential error estimates for the metric problem 
obtained in the previous section with a comparison argument to obtain estimates on the difference 
between —5v s (y,ui ;p) and H(p). The comparison argument, which was introduced in [4] to prove 
homogenization, yields an estimate from below for 5v s + H(p) for all p € M. d and from above only 
for p's away from the flat spot. 

Proof of Theorem [2](i). We begin with the demonstration of (|2.15p . Fix < 5 < X < 1 and 

uj G n such that -5v s (0, uj ;p) > H(p) + A. By (J2J2) and (JO}, for < c < 1 sufficiently small, 

w(y) :=v 5 (y,uj;p)-v s (0,uj;p) + cX(l + \y\ 2 p - cX 

satisfies 

H(p + Dw,y,uj) > -5v 5 (y,uj; p) - -X in R d . 

According to f)5.2|) . there exists y\ G R d such that \yi\ < C/X5 and 

w(yi) = inf w(y) < w(0) = 0. (5.4) 
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Define U := {y 6 M d : w(y) < A/4<5} and notice that, for every y € U, 

-6v s (y,oo;p) > -5v s (0,u> ;p) - ^A > H(p) + jX. 

In particular, 

H(p + Dw,y,ui) >H(p) + -A in U. 

i 

Next we denote w(y) := w(y) — w(yi) + c\ (l + \y — yi| 2 ) 2 — cA and observe that, by \Dw\ < K p + 1 
and (|2.5p . if c > is small enough, then 

H(p + Dw,y,uj)>H(p) + \\ in U. (5.5) 

Also notice that w > w, w(yi) = and w{y) > cA (l + \y — yi\ 2 ) 2 — cA. Therefore, if we define 
V := {y G R d : w{y) < A/45}, then, for some < R < C/5, we have 

VCUnB( yi ,R). (5.6) 

Let fj, := H(p) + A/4 and m(y) := — m M (yi, y, oj)—p-{y — y\) and observe that, according to (|A.14|) . 

H(p + Dm,y,oo) <H(p) + ^X in R d . (5.7) 
The comparison principle (Proposition 13.41) yields 

max (-m M (yi, y, w) +p • (yi - y)) = max (m - w) + A > m(y x ) - w{yi) + — = — . (5.8) 
ySidV aV 40 40 40 

Since > H(p) we see from (IA.38D that p 6 dfn^(0), that is, for every z 6 M d , 

^n(z) > p • z. (5-9) 

Using (pTgjh and the fact that dU C B( yi ,R), we obtain 

inf (m M (yi,y,w) -m M (yi -y)) < -A 

yeB( yi ,R) 4d 
In view of (jA.lip and |yi| < C/X5, we deduce that, for some c > small enough, there exists 

?C/A<5> 



z £ K := {cXk/5 : k e Z d } n B c 



such that 

inf (m fi (z, y ,co)-m fl (z-y)) <-—. 

y£B(z,R) OO 

In summary, we have shown that, for some fixed R < C/5, 

LeH:-fo 5 (0,w iP )>%)+Ak M Lefi: inf (m^z, y, u) - m^(z - y)) < -A j . 
I. J ^j^. ^ y<=B(z,R) 00 J 

Hence, using (HI]), 5 < A, that |Jf| < CA~ M , p > A/4 and (|435l) . we obtain 

A " 



-5v d (0,-;p) >H{p) + X 



< CX~ 2d l 



inf (m^(y,0,-) - m jU (y)) < 

<Cr^e xp (-^W«e x p(-^ 



This completes the proof of (|2.15p . □ 
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Proof of Theorem [2](ii). The proof of (|2.17p is similar to that of (|2.15p . Since it is a bit more 
complicated, we give complete details. Fix 00 E f2 for which — 5v s (0,uj ;p) > H(p) + A. According 
to (I2.5p . (|5.3p and Lemma [331 if c > is chosen sufficiently small, then 

w{y) := v s (y,uj;p) - v s (0,uj;p) - c\\y\ 

satisfies 

H(p + Dw,y,u) < -6v s (y,uj;p) + -A in R d . 

Moreover, according to (15. 2p . there exists y 2 E K rf such that 1 2/2 1 — C/XS and 

w (y 2 ) = S up > io(0) = 0. (5.10) 



Define U := {y E M rf : > —A/45} and observe that, for every y E U , 

-5v 5 (y,uj;p)<-5v 5 (0,uj;p) + ^X<H(p)-^X 

and, therefore, 

H(p + Dw,y,uj) <H(p) - -A in {7. 

By |Du>| < C and (I2.5p . there exists c > is small enough that w(y) := u>(2/) — w{y2) — cX\y — y 2 \ 
satisfies 

H(p + Dw,y,uj) <H~(p) - -X mU. 

Let V := {y € M d : u5(y) > —A/45} and note that, in light of the fact that w < w, 10(2/2) = and 
w{y) < — cX\y — 2/2 1) there exists D < R < C/d such that 

V cUnB(y 2 ,R). (5.11) 

For notational convenience, denote [x := H{p). Since H(p) > 0, according to (|A.39P there exists 
e E R d with |e| = 1 such that p E 9m^(e), that is, for every £ E IR rf , 

m^z) — p ■ z > = m^(e) — p • e. (5-12) 

Define 777,(2/) := m^(y, y 2 — Re, uj) — p ■ y and observe that 

H(p + Dm, y, uj) = H(p) in R d \ {y 2 - Re} D V. 

An application of the comparison principle yields 

w(y 2 ) - m(y 2 ) < sup (w(y) - fh{y)) = max (w(y) - rh(y)) . 

Using (|5.1U|) and that w(y 2 ) = and w = —A/45 on dV, we obtain 

max (m^(y 2 , 2/2 - Re, uj) - m^y, y 2 -Re,u))-p - (y 2 - y)) > max (777(2/2) - 777(7/)) > 

y£dV y&dV 45 

It follows that either 

m^(y 2 ,y 2 - Re,uj) -fn^(Re) = m tl (y 2 ,y 2 -Re,uj)-p- (Re) > ^ (5.13) 

or else 

max (-m^(y, y 2 - Re, uj) + m M (y - 2/2 + -Re)) 
y€av 

> max (-m M (y, 2/2 - Re,uj) + p ■ [y - y 2 + Re)) > (5.14) 
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Using (jA.lip . ()5.1ip and that 1 2/2 1 < C/\8, we see that there exists 

z e K := {cXk/5 : k G Z d } n B c/X5 , 

such that either 

m^(z, z - Re, uS) - m^(Re) > 

or else 



sup (-m M (y, z - Re, u) + m^(y - z + Re)) > —. 

y&B(z,R) °° 
In summary, we have shown that 

[uj^VL : -Sv 6 (0,oj;p)>H(p) + x\ C [j (E^z) U E 2 (z)) 



(5.15) 
(5.16) 

(5.17) 



where E\(z) is the event that (|5.15p holds and E 2 (z) is the event that (|5.16p holds, and R < C/5 
is a fixed constant. Using (|A.4j) and \K | < CX~ 2d < C5~ 2d , we obtain 



-5v s (0,-;p) >H{p) + X <C5- 2d {P[E 1 (Re)]+¥[E 2 (Re)]). 
Applying (I4.55p . we find that 



W[E 2 (Re) 



inf (m M (y,0,u;) -m^y)) < 
y£B R 



A 



< Ci? d exp 



C5 2 R 



< C5~ d exp 



(5.18) 
~~C5 



and, if A satisfies (|2.16p for sufficiently large C3 > 0, then we may apply (I2.14p to get 



•[Ei(Re) 



m^(Re, 0, •) - m^(Re) > 



< exp 



pX 2 
C5 2 R 



< exp 



pX 2 
~C6 



Combining the last two sets of inequalities with (I5.18P yields (I2.17p . completing the proof. □ 

A covering argument now yields explicit error estimates for (|5,ip in balls of radius O (8~ 1 )- 

Lemma 5.1. For each K > 0, there exist constants C\q,c 2 > such that, for each p E Bk, R> 
and < 5 < c 2 , 



sup -8v (y, ■;p)> H(p) + Ci 6 | log 5| 3 <J; 



< R d 5 2 , 



(5.19) 



and, if H(p) > 0, then 

inf -5v s (y,-;p) <H(p)-C 16 (~H(p)-%8^ +H{p)~ l 5^) (l + |log<5| + \\og H{p)\)* 



< R d 5 2 . (5.20) 



Proof. Since both of the estimates are obtained from the first two statements of Theorem [2] using 
a similar argument, we prove only (|5.19p . To do so, we apply (|2.15p with A := A\ log 5\ 3^3 , for 
A > chosen sufficiently large, and use a simple covering argument. Notice that if < 5 < \ is 
sufficiently small, depending on A, then 5 < X. There exist points y\,...,yN € B R /g such that 
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N < CR d \- d and the balls B(y j ,X/2SK p ) cover B R/S . According to (lETo]) . (O]) and ([2Tl5j) . with 
C3 > as in (j2. 15j) . we have 



sup -8v 5 (y, ■ ;p) > H{p) + A 



< 



max -5v 5 (yj, ■;p)> H(p) + - 

l<J<iV Z 



< CJVH 



-6v s (0,-;p) >H(p) + 



A 



8C 3 <5 
A 3 |log<5| 
8C7 3 

We therefore obtain ([5TT9J) if we choose ,4 > so that A 3 > 8C 3 (4d + 2). 



< CR d \- d 5- 3d exp 

< C^-^exp 



□ 



We next apply Lemma [57T1 along a certain subsequence 5 n — > to prove, with the help of (IB.lOj) 
and the Borel-Cantelli lemma, the first two statements of Theorem [3j 



Proof of Theorem [3](i) and Theorem [3](ii). The arguments for the two statements are almost 
identical, so we prove only (123UD . Let R > 0, 5 n = n' 1 and apply (ETT9D to obtain 



E 

n=l 



sup -5 n v 5n (y, ■;p)>H(p) + C 16 | log S n | s 5, 

V^ B R/S n 



n>l/C2 



OO. 



n=l 



By the Borel-Cantelli lemma, we deduce that, for each R > 0, 



lim sup sup 



-S n v Sn (y,-;p) - H(p) 



< 1 



y^i^ Ci 6 1 log cykl 



Intersecting these events for each R = 1, 2, 3, . . ., we find an event of full probability such that, 
for every R > and ui £ Q\, 



lim sup sup 

n— >oo 



-5 n v s "(y,-;p) -H(p) 
1 

r3 



< 1. 



y^ B ^n Ciellog^li^ 
Notice that 5 n+ i/5 n = 1 — 5 n +i and so, according to (jB.lOp . for any S n+ % < 77 < S n and oj G Q, 



sup 



5 n v Sn {y,uj;p)-'qv r, {y : ixi;p) 



< C5 n+ i < CS n . 



Hence for every w S Oi we have 



-6v 5 (y,u;p) - H(p) -r)v T >(y,u ;p) - H(jp) 
urn sup sup j — j = lim sup sup sup j — j 

<5^0 yeB R/s Cie\log5\3Ss n^oo 8 n <r)<8 n ^ 1 yeB R/v Cl6 | log 77 1 3 77 3 



< lim sup sup 

n->-oo y€B 2R/5rl 



-5 n v^(y,u;p)-H(p)+C5 ri 
1 1 

C 16 \\og6 n \'36% 



< 1. □ 
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5.2. The sub-ballistic regime. We show that the behavior of 5v s (0, • ;p) for p's on the flat spot 
{H = 0} is determined by the distribution of H(0,0, •) near its maximum and that with further 
(quite reasonable) assumptions on this distribution, we obtain exponential error estimates and an 
algebraic rate of convergence for — 8v s (0,u ;p) to H(p). 

We begin with the simple observation, which is probably well-known and essentially taken 
from [3], that —5v s (-,u ;p) is controlled pointwise from below by H(0, -,u>). 

Lemma 5.2. For every p E M. d and uj £ U, 

— 5v s (0, uj ; p) > sup sup H(0,y,u>) — K P R5 , 



(5.21) 



R>0 \y£B R 

with the constant K p > defined in (|B.6|) . 

Proof. Fix p,y G R d and u> G fi. Due to (H3D and (^9]) . w 5 (-,o; ;p) satisfies 

Sv s + H(0,y,u) < 5v 6 + H(p + Dv s ,y,u) = in E d . (5.22) 

While this holds in the viscosity sense, there are no derivatives in the expression on the left of (|5.22p . 
and we deduce that 

-Sv\y,u]p)>H(0,y,u). (5.23) 
Combining QEggJ) and J53D yields (I5^TD . □ 

For p's on the flat spot, the rate of convergence for — 5v s (0,uj ;p) to = H(p) given by (|5.2ip 
is essentially optimal (see Lemma 15.51 below). We next exhibit exponential error estimates for 
—5v s (0, ■ ;p) under the additional hypothesis (12. lip . 

Proposition 5.3. Assume (|2.1ip . There exist Cyj, C3 > such that, for all 8 > and < A < C3, 

X d+e ■ 



-Sv d (0,-;p) < -A 



< exp 



C 17 8 d 



(5.24) 



Proof. Fix p£R d and apply (fOTT) with R := X/(2K p S) to discover that 



-<5u*(0,-;p) < -A 



< 



sup H(0,y,-) < K p R5-X 

y^B R 



sup H(0,y r ) < ~ 
y&B R 2 



Select yi, . . . , yjv € -Br such that iV > cR d and, for every i 7^ j, |yj — > D. Then (|2.2|) and (|2.3p 
yield 



sup H(0,y,-) < ~ 
y&B R 2 



< 



sup H(0,yj, •) < ~ 
l<j<N * 



n > 

l<j<AT 



< 



A 



Using (|2.1ip . we find 

H r\H(o, yj ,-) 



< 



#(o,o,o < 



A? 



< [1 - c\ 

< exp f-cNX 6 } < exp f-cA 



d+0 fi-d 



Combining the lines above yields (|5.24p . 



□ 
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We now obtain, under assumption (|2.1ip . exponential error estimates and a rate of convergence, 
from below, for —5v s for all p £ M. d . First, we combine (|2.17p and (|5.24p to obtain error estimates 
independent of H(p). 

Proof of Theorem [2](iii). Observe that if we take A > sufficiently large and a > sufficiently 
small, then, for any < 5 < a and A > satisfying (j2. 18[) . we have 



A > A6*\\og5\* > sup (2a A 4C 2 (a~2§2 + a ' L 5^ (1 + |log<5| + |log<r|)2 . (5.25) 

cr>0 ^ ^ / / 

Thus (I2.18P ensures that ()2.17p and (I5.24p overlap in an appropriate way to yield the theorem. 
Indeed, for < A < C3, we may apply (|5.24p in case A > 2H(p), while in the case that A < 2H(p), 
we may apply (|2.17p . since (|5.25p ensures that A < 2H(p) implies (|2.16p . □ 

Arguing in a similar way as in the proof of Lemma 15.11 we obtain the following result as an 
application (|2.19p . The details are left to the reader. 

Lemma 5.4. Assume (|2.1ip and let a and ft be defined as in (I2.22p . There exist constants C\g, C4 > 
such that, for each R > and < 5 < C4, 



inf -5v 5 (y,-;p) < H (p) - C 18 \log Sf 6° 
Using (|5.26p we complete the proof of Theorem [3j 



< R d 5 2 . (5.26) 



Proof of Theorem [3](iii). The statement follows from (|5,26p by an argument very similar to the 
proof of Theorem [3ji) given above. □ 

5.3. The rate may be arbitrarily slow on the flat spot. As explained in Appendix the 
vector p = belongs to the flat spot: H(0) = 0. We show here that (|2.1ip is necessary for the 
existence of an algebraic rate of convergence like (|5.26p for the limit (|5.ip at p = 0. Furthermore, 
without some assumption on the distribution of the random variable H(0, 0, •) near its maximum, 
there is no restraint on how slowly the limit (|5.ip may converge for p = 0. 

We begin by exhibiting an upper bound for — 5v s (0, uj ; 0) to match (|5.2ip . which shows that, for 
p's on the flat spot, the rate given in Lemma 15.21 is essentially optimal. 

Lemma 5.5. There exists a constant C19 > 1 such that, for every uj € Q and R, 8 > 0, 

r) Ft 

- 5v 5 (0, uj ; 0) < {-SR) V - sup H(0, y, uj). (5.27) 

C-19 + OH y£B R 

Proof. Fix uj € f2, R, 5 > and define, for every < a < 1, 

h R := - sup H(0,y,uj) and h a ^ R :=- sup H(p,y,uj), 

y£B R ( P ,y)eB a xB R 

where we set h a ,R '■= h R in the case that a = 0. Observe that h R > and, due to (|2.5p . for any 
< a < 1, 

\hn - h aiR \ < Ca. (5.28) 

Next fix < a < 1 to be selected below, define u(x) := a(R — \x\)+ and note that u € L, \Du\ < a 
in R d and u = in R d \ B R . We conclude that 

8u + H(Du, y, uj) < (SaR - h a , R ) V in R d . 
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Choose < a < 1 to be the largest number for which a5R < h aj R A 5R. The comparison principle 
(Proposition EH) yields that 5aR = 6u(Q) < 5v s (0,w ;0). According to ([OS]) . 



a > 



C + 5R 



A 1. 



which concludes the proof. 



□ 



The previous lemma states that the rate — 5v s (0, u ; 0) converges to is controlled from below, 
up to a factor of 2, by the maximum of H(0,-,ui) in the ball B Cig /$. Using the independence 
assumption and an easy covering argument, we relate the latter to the distribution of H(0, 0, uj) 
near its maximum to recover the following estimate. 

Lemma 5.6. There exists a constant C20 > such that, for every 5 > and < A < 1, 



logP \-6v s (p,-;0) < -A > C 20 S- d log ^1 - C 20 A~ d P[fl'(0,0,-) > -4AM , (5.29) 

where the inequality is vacuous if the argument in the logarithm on the right side is negative. 

Proof. Fix < A < 1. There exists a finite collection of points {yij : 1 < i < N,l < j < M} C Br 
such that < CR d , M < C\~ d , Br is covered by the balls B(yij,cX) and \yij — ykj\ > D if i ^ k. 
According to (^2]) . (|2~3|) and ([23]), we have 



sup H(0,y,-) < -2A 
y£B R 



> 



sup sup H(0,yij,-) < — 4A 

l<i<iV 1<7<M 



sup H{0, yij ,-) < -4A 

i<i<Af 



N 



> (1 - MP [#(0,0,0 > -4A]) 



N 



Setting R := Cig/5, with C19 > 1 defined in Lemma [531 applying (|5.27p and taking the logarithm 
of the resulting expression yields (|5.29[) . □ 



We next show that the assumption (|2.1ip is essentially necessarily for (|5.ip to have an algebraic 
rate of convergence at p = 0. 

Proposition 5.7. Assume, contrary to (|2.1ip . that there exists C, c > and 9 > d such that, for 
every < A < c, 

P [#(0,0,0 > "A] < C\ e . 
Then there exists C5 > such that, for every < 5 < C5, 

-6v s (0,-;0) < -<5 7 1 > c 5 /or j:=d/(0-d). 



(5.30) 
(5.31) 



Moreover, for uj belonging to an event of full probability, we have, for every rj > 7, 

<fo*(0,w;0) 



lim inf ■ 

5-s-O 



+00. 



(5.32) 



Proof. Let A > 0. Using the elementary inequality — 2t < log(l — t) for < t < i, we apply (|5.29|) 
and (|5.30p to obtain, for sufficiently small 5, A > 0, 



log I 



-<Ju tf (0,-;0) < -A 



> C5- d \oz ( 1 - CX e - d ) > -C5- d \ - d . 



Setting A := 5^ yields (I5.3ip for sufficiently small 5 > 0, while on the other hand setting A = 5 a 
for a > f3, applying the Borel-Cantelli lemma along a subsequence and arguing as in the proof of 
Theorem [3] yields (1021) . □ 
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It is clear from (|5.29p that, if f[H(0, 0, •) > —A] decays to very quickly as A — > 0, then 
-Sv 5 (0,u;0) will converge to zero very slowly. We conclude this section by exhibiting an example 
demonstrating an arbitrarily slow rate. 

Example 5.8. Let p : [0, 1] — > [0, oo) be a given increasing continuous function with p(0) = 0. We 
will construct a Hamiltonian H, satisfying (|2.1U|) . such that 

lim inf ^ V ^ ' ^ > 1 almost surely. (5.33) 
<5^o p{5) ~ y ' 

According to ()5.29j) . it suffices to exhibit an H for which 

C 20 log (l - C 2 op(5)- d P [H(0, 0, •) > -4p(<5)]) > -5 d+2 . (5.34) 

Indeed, (IQ9|) and (Q4l yield 

P \Sv s (0, • ; 0) > p(S)\ > exp(-C5 2 ) >1-C5 2 

for 5 > sufficiently small. Then the Borel-Cantelli lemma and an argument similar to the proof 
of Theorem [3j using (|B.10p , give (|5.33j) . 

It is easy to check that for (|5.34j) . it suffices to have 

P [H(0, 0, •) > -a) < ca d (p-\a)) d+2 =: p(a), (5.35) 

and so we need to construct a Hamiltonian with a very thin distribution near its maximum. This 
is quite simple. We may take, for example, 

H(p,y,cj) : =lH 2 -0(F( y , w )) 

where V is a Poissonian potential (see for example [33]) and (p : [0, oo) — > (0, oo) is a continuous, 
decreasing function such that (p(t) decays very slowly to as t — > oo (the precise rate of decay 
required can be explicitly calculated in terms of jo). We leave it to the interested reader to fill in 
the details. 

5.4. Uniform error estimates for the approximate cell problem. The proofs of Theorems H] 
and El given in the next section, depend on the following extensions of Theorems [2] and [3] which 
hold uniformly for bounded \p\ and for y in balls of radius 0(<5 -Ar ), for any N > 1. We omit 
the arguments, since the error estimates follow easily from an application of Theorem [2] combined 
with (|B.5j) and a routine covering argument, and then the convergence rates from the latter using 
the nearly same argument as in the proof of Theorem [3l 

Proposition 5.9. For each K > 0, there exist constants C2i,cq > such that, for each R > 0, the 
following hold: 

(i) For every < A < 1 and < 5 < cqX, 
sup (-6v s (y, ■ -p) - H{p) \ > A < C 21 R d 5~ 5d e W (-^) • (5.36) 

_(y,p)€B R xB K V ' J V <-21<V 

(ii) // (I2.11|) holds and A, 5 > satisfy < 5 < c§\ and 

A > C 2 i<5s| log<5|3, (5.37) 
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then 



inf 

(y,p)£B R xB K 



-5v 5 (y,-;p)-H(p)) < -A 



< C 2 iR d 6- 5d exp I --- I '— A 



1 /A 3 A dH 



C, 



21 



5<< 



. (5.38) 



Proposition 5.10. There exists an event SI3 S J of full probability such that, for every u G Q3, 
the following hold: 



(i) For every R,K>0 and N > 1, 



-dv 5 (y,uj;p) - H{p) 



limsup sup sup 

<5^o yeB R/gNP eB K 5^ | log c5| 



< 00. 



(5.39) 



(ii) // (|2.1ip holds and we define a and f3 as in (|2,22p . then, for every R, K > and N > 1, 



lim inf inf sup 

5^0 yeB R/sN p&Bk 



-5v s (y,uj;p) - H(p) 
5 a \\og5\P 



> -co. 



(5.40) 



6. Error estimates for homogenization 



We now present the proofs of Theorems H] and [5j The main step is to precisely quantify how 
the 5v Sl s control the u e 's, so that we may apply the results of the previous section to obtain error 
estimates and a rate of convergence for the latter. 

For each e, T > 0, let u e = u £ (-, -,u>),u = u(-, •) <E BUC(lR d x [0,T]) be the unique solutions of 

u £ t +H (du £ ,-,uj) =0 in R d x (0,oo), f u t + H (Du) = in R d x (0, 00), 

V e J and < 

u £ (-,0,oj) = n , [u(-,0)=n , 

where u £ C ' 1 ^) is the given initial data with 1 1 ^0 1 1 c ' 1 (IR. 1 *) — 

It turns out (see, e.g., [6]) that there exists a constant L > 0, depending on K and the assumptions 
for H, such that, for all e > 0, x, y £ R d , s,t>0 and well, 

\u £ (x,t,uj) -u £ (y,s,Lo)\ < L(\x-y\ + (6.1) 

and 

\u{x,t)-u{y,s)\ <L{\x-y\ + \s-t\). (6.2) 

These estimates are derived principally from the coercivity of H. Recall that, due to (|5.2p 
and (|B.12p . the effective Hamiltonian shares the same rate of coercivity assumed in ()2.6|) . It 
also follows easily from this that for each e > 0, x £ M d , < t < T and oj G fi, 

|«(x,t)| + |u £ (x,t,w)| <K + LT < C(l + T). (6.3) 

The important link between the 5v s, s and the -u e 's is summarized in the following lemma. Then 
Theorems |4] and [5] follow relatively easily from it and Propositions 15.91 and 15.101 The basic idea 
is that the event that \u £ (x,t,u)) — u(x,t)\ is large should only be observed if \5v s + H{p)\ is also 
large. The proof, which is rather technical and lengthy, follows along the lines of [8] with necessary 
modifications to deal with to the lack of uniform estimates on the difference between — Sv s and 
H(p). It essentially consists of quantifying the perturbed test function method [13] to argue that, 



ERROR ESTIMATES AND CONVERGENCE RATES FOR STOCHASTIC HOMOGENIZATION 35 

if —5v s is close to H(p), then, up to an appropriate error, it properly captures the oscillations of 
u e , that is, 

u(x, t) « u £ (x, t, uj) — ev 5 ^— , uj ; Du(x, t)j , e <C 5. (6.4) 

Rather than apply the comparison principle, we must use the proof of it. This is due to the fact 
that, since u is not in general C 1 , we cannot insert p = Du(x,t) into v s (x,uj ;p). There are other 
technical difficulties (the presence of three nonsmooth functions, the fact that v s is not smooth in 
p) which we handle by the standard viscosity theoretic technique of doubling (or rather tripling) 
the variables. 

Lemma 6.1. For each K > 0, there exist C22,L > such that, for every X,T,e,5 > satisfying 

A < 1, T > 1, e < XT and < 5 (6.5) 

TX A 

and each initial datum uq € C 0,1 (]R rf ) with \\uo\\ c o,imd\ < K, we have 



uj G : inf inf (u e (x,t,uj) — u(xA)) < — C22AT 
xeB T 0<t<T 



C sup sup [ -5v (y , uj ; p) - H (p) ) > A f , (6.6) 

I |s/|<C 22 5AT 2 /e 2 |p|<£ 



and 



uj € VL : sup sup (u £ (x,t,uj) — u(x,t)) > C22AT 
x&B T 0<t<T 



C iuen : inf .inf (-5v a (y,u> ;p) - H(p) < -A . (6.7) 

L \y\<C220\T z /e 2 \p\<L \ /J 

Proof. Since the arguments for (16.6h and (16. 7h are nearly identical, we prove only (|6.7p . 
Define C : K d -> M d by 

L A |x| 

and notice that C(x) has the same direction as x and 

\({x)\ = LA\x\ and |C(a?) - C(y)| < \x - y\. (6.8) 
Fix T, A, e, 5 > satisfying 

A < 1, T > 1, e < AT and — r < 5 (6.9) 

TX Z 

where the constant A > 2 is selected below. Also fix parameters a, 7 > to be chosen below and 
consider the auxiliary function : M d x M d x [0, T] x [0, T] x — > R given by 

$(»,y,t 3 s,w) := u e (x,t,w) - u(y, s) - ev 5 fp w ;C ~ ^~' x ~ y ' 2 ~ ^~ (* ~ s ) 2 

-As-7(l + |x| 2 )^+7. (6.10) 
Using ||521> . (1531) and (|5TH]) . we have, for each (x,y,t,s,w) € R d x M d x [0,T] x [0,T] x O, 

\§(x,y,t,s,u)\ < C{1 + T) + Ced' 1 - -L| x - y| 2 - ^-{t - s) 2 - As - 7 (l + |x| 2 )^ + 7. (6.11) 

la le 
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It follows that, for each u> € fi, the function 3>(-,w) attains its global maximum at some point of 
R d x R d x [0,T] x [0,T]. Set 

M(w) := max 3>(-,w) 

M d xR d x[0,T]x[0,T] 

and denote by E £ 3" the event that the global maximum of <£(-,u;) is attained by some point 
(x, s) with either t = or s = 0; that is, 

£:= w£Sl:M(w)= sup <l>(x, y, 0, s,oj) V <&(x,y,t, 0,lu) > 

[ (z,j/,t,s)eR d xR d x[0,T]x[0,T] J 

Note that, for notational simplicity, we omit the dependence of E on a and e. 

The lemma follows from the fact that, under an appropriate choice of the parameters, 

sup sup sup (u £ (x, t, u) — u(x, t)) < CAT. (6-12) 

ueExeB T o<t<T 

and, for some R < C5XT 2 /e 2 , 

0\£cLg(] : inf ( -5v S (y,oj ;p) -H(p)) < --] . (6.13) 
I (y,p)€B R xB L \ / 2 J 

Proving (16. 12ft is relatively easy, while (I6.13P involves a more involved comparison argument. 

Before we begin the proof of (16.130 . we record two elementary estimates that necessarily hold 
at any global maximum point of and for any uj £ f2. For the moment, we fix oj € O and a 

point (x ,y ,t ,s ) Gi^xl^x [0,T] x [0,T] satisfying 

<S>(x ,y ,t ,s ,Lu) = M{oj). (6.14) 

The inequality $(xq, yo, to, sq, oj) > <J>(0, 0, 0, 0, w) yields, in light of (|6.1ip . that 

1 Ce 
7ko| + tt^o - so) 2 < C(l + T) + — < CT. (6.15) 
2e o 

If so 7^ 0, then by (|6.2j) and that s i— > u(y, s) + (s — t) 2 /2e has a local minimum at s = sq, we deduce 
that 

[s — *o] < (L + X)e < (L + l)e. (6.16) 

The inequality (16. 16H is also satisfied for a similar reason if to 7^ 0, and trivially if so = to = 0. We 
also claim that 

\x - yo I < La. (6.17) 
If not, then y i— > C(( x o — y)/ot) is constant in a neighborhood of yo and we obtain from (|6.14p that 

y h-> u(y, so) H |xo — y\ 2 has a local minimum at y = yo- 

2a 

Due to (|6.2p . we conclude that |xq — yo I < Ta, which is a contradiction. In particular, 

c (xo^yo\ = xo-yo^ 
\ a J a 

We now commence with the proof of ()6.13|) . following the classical argument of the proof of the 
comparison principle for viscosity solutions and [8j . Fix uj G Sl\E and select a point (xq, yo,to, so) € 
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R d x R d x (0,T] x (0,T] such that (foTTil) holds. The first step is to fix (x,t) = {x ,t ), allow (y,s) 
to vary, and to use the equation for u. The goal is to derive ()6.23p . below. From (|6.14j) . we see that 

5 ( x . > f x ~ V\\ , 1 | i2 , 1 n \2 



(y,s)^u(y,s) + ev d (->; C (^n + — |x - y| 2 + — (t - s) 2 + As 

has a local minimum at (y, s) = (yo> sq). (6.19) 

According to ([EH and (fBTgj) . 



5 / a; / x - y \ \ 5 / x / x - yo 
— ,w;C - £v — > W ;C 



e \ a I I \ e \ a 



Ce 
< — 
- 5 



x Q -y \ _ / x - y 
a I \ a 



< Ce \y-y°V ( 6 . 2 o) 

oa 



Using (|6.19p , (|6.20p , the fact that equality holds in ()6.20p at y = y and by enlarging C > slightly, 
we obtain that 

(y,s) (->■ u(y,s) + —\x - y| 2 + — (t - s) 2 + As + C^-\y - y Q \ 
la 2e oa 

has a strict local minimum at (y, s) = (yo, sq). (6.21) 

It follows that, for all sufficiently small /3 > 0, there exist (yg, sp) S M rf x [0, T] such that (yp, sp) — > 
(yo, so) as /3 — > and 

(y, s) h+ n(y, s ) + ±-\x - y| 2 + —(to - s) 2 + As + C^- (/? + |y - y | 2 )^ 

has a local minimum at (y,s) = (yp,sp). (6.22) 

Using the equation for u, we obtain 

-\+~(h-sp) + H(^^-Qp\>0, 

_ i 

where := C^j (/3 + [t/^a — yo| 2 ) 2 (y/3 — yo)- Since \Qp\ < Ce/5a, the Lipschitz continuity of iJ 
yields 

and, after letting /3 — >■ 0, we find 

- A + i(to - so) + 5 f > -C-I-. (6.23) 

e \ a J oa 

We next fix (y, s) = (yo, so) and let (x, t) vary, in order to use the equations for u £ and v s . The 
intermediate goal is to prove (|6.28|) . below, to compliment (|6.23[) . This leads to some analysis that 
is a little more complicated that what we have just performed above to produce (I6.23j) . due to the 
fact that we have two (in general, nonsmooth) functions u £ and v s which depend on the variable 
x. This sort of technical difficulty is typically handled using an "iterated" perturbed test function 
argument (an idea introduced in [13]). We need to quantify this idea, and so, following [8j, we 
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fix another parameter a > (which will be sent to zero shortly) and introduce a second auxiliary 
function * : R d x R d x [0, T] -> R defined by 

V(x,z,t) := u £ (x,t,u) - ev 5 (^, u -x(^^X\ - ±-\x - y \ 2 - ^(t - s Q ) 2 

-7(l + |x[ 2 )^ -— \z - x\ 2 - j (l + \x - x \ 2 )K (6.24) 

The last term in (|6.24p provides some strictness and therefore, by (|6.14p . there exist points 
(x CT , z a , t a ) G M. d x M. d x [0,T] such that (x a , z a ,t a ) — >• (xo,xo,io) as a — > and 

^(x^,^,^) = sup 

R d xR d x[0,T] 

Freezing z = z a and letting (x, i) vary, we have 

(x,t) ^ u £ (x,t) - — \x - y \ 2 - — (t - so) 2 - 7 (1 + k| 2 )^ - ttW ~ x\ 2 - j (l + \x - x | 2 )^ 
2a 2e 2<r 

has a local maximum at (x,t) = (x a ,t a ). 

It follows from the equation for u £ that 

-(*<r - S ) + Hi 1 \~ P(j , ,U)\ < 0, 

e \ a <t e J 

where P a := 7 (l + \x a \ 2 ) 5 x CT + 7 (l + \x a — xo| 2 ) 5 (x CT — xq). Since \P a \ < Cj, we use (|2.5p to 
obtain 

i(t (T - So ) + ^f^^ + ^^,-,^ <C 7 . (6.25) 
On the other hand, freezing (x,i) = (x ff ,f ff ) and letting z vary, we get 

2 1 — y ev s ( — , uj ; C f — ] I H |z — xJ 2 has a local minimum at z = z^. 

V e V a ) ) 2d 



According to ([6.8 



■ z ~ 2/0 \ _ . / z CT - 2/o 



<— z-zJ. (6.26) 



a J \ a 

with equality at z = z a . Hence by using (|B.9P again and enlarging C slightly, we obtain 

. 5 f z j- ( Zv-Vo \\ . 1 1 ,2 , n £ I I 

z 1— > £D l-,w;CI — - — I I + — |z - x °\ +C ~fa) z ~ ZcT \ 

has a strict local minimum at z = z a . 

Thus we can find points z a>K — > z a as k — > 0, such that 



e \ a J J 2a 5a 

has a local minimum at z = z CTift . 

Using the equation for v s , we discover that 

6v* w ; C (^) ) + H (C (^) + - P„„, ^, a,) > 0, 
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_ l 

where P a<K := Cj^ [n + \z a , K - z a \ 2 ) 2 (z a>K - z a ). Since \P a , K \ < Ce/5a, we use ([23]) to get 

\ e \ a J J \ \ a J a e J oa 

Letting k — > yields 

A, (* „ if (*ZJ»)) +H ( ( (iZ") + * „) > - C A. ( , 27) 

Comparing (j6.25j) and (j6.27j) . using (j2.5j) . (jB.9[) and (16.26j) and then sending a — > yields 

^ ( ^,0, ; X -^] > -(to - -o) - Cf - C7 7 . (6.28) 
\ s a J e da 

We now combine (|6.23j) and (16. 28ft to deduce, upon choosing 7 := e/(8a), that 

-ft,' f^j^ZW^Hf^^Ucf -A. 
\ e a J \ a J da 

Select a := AT to deduce that, for every 5 > Ae/TX 2 = Ae/(aX), we have Ae/5a < A, and hence 

-5v s (^■^^] <h(^^]-Ix 



a J \ a J 2 
provided that A > 2 is chosen sufficiently large. In light of (16. 15f) and (16. 17H . we have shown that 

n\Ec\u>en: inf ini (-6v 6 (y,u ;p) -H(p)) < -~x\ , (6.29) 
[ y£B R p£B L \ ) 2 J 

where 

bo I CT C5XT 2 

R-= — < — = — — , 

e 76 e z 

and we have proved (|6.13j) . 

We turn our attention to (f6Tl2jh Note that 5 > Ae/TX 2 > Ae/TX and 27 < A since A > 2. 
Therefore, for any w £ fi, 1 £ Bt and < f < T, 

u e (x,t,uj) - u(x,t) = $>(x,x,t,t,U)) + ev 5 (x,oj; 0) + At + 7(1 + \x\ 2 )^ - 7 

< M(u) + CefT 1 + AT + 7 T < M(u) + 3AT. (6.30) 

Fix to 6 £ and let (x ,yo,t , s ) £ M. d x R d x [0,T] x [0,T] be a point which satisfies (flTlIj) and 
such that either to = or so = 0. If to = 0, then using (|6,16p . (|6.17p . (|6.9p and the assumption 
that < e < XT, 

M(u>) = ®(xo,yo,to,s ,uj) < u (x ) - u(y ,s ) - ev S ( — ,u ; — — — 

\ e a 

< L\x - yo \ + Ls + Ce5~ l < L 2 a + C(L + l)e + Ce5~ l < CXT. 

If so = 0, we also get the bound M(oo) < CXT via a similar observation. Combining this with (|6.30|) . 
we obtain 

sup sup (u £ (x,t,uj) — u(x,t)) < CXT. 

x&B T 0<t<T 

This completes the proof of (|6.12p and thus that of (|6.7j) . □ 
We are now ready to complete the proofs of our main results. 
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Proof of Theorem [4j The theorem is easily obtained from the combination of Lemma 16.11 and 
Proposition 15.91 For example, to prove (i), we take 5 := C 22 e/TX 2 and check that A > Ces for 
large enough C > and T > 1 imply that < 5 < cqX. Thus we may apply (|5.36p to (|6.6[) . To 
prove (ii) , we define 5 in the same way and check that (|2.27[) implies (|5,37p so that the combination 



of (I5.38P and (|6.7p yields the estimate. □ 

Proof of Theorem [5]. The result follows from a combination of Lemma l6.1l and Proposition 15.10] 
We first prove (i). Define, for each e > 0, 

A(e) := Ae^\logefi and 5(e) := C 22 e/TX(e) 2 , 

where A > 1 will be selected below. Notice that, for large A, 

A(e) > cA 2 5(e)\ |log<5(e)|3 > A5(e)\ |log$(e)|s . 

According to (|6.6p . 



P| |J jwGft: inf inf (u £ (x,t,uj) -u(x,t)) < -C 22 TA(e) 1 

V>0 0<e<ri { xeB T 0<t<T J 

C n M L G O : sup sup (-5(e)v 5 ^ (y, w ; p) - #(p)) > A(e) I 

r?>0 0< £ <^ I |j/|<C 2 25( £ )A( £ )T2/ £ 2 bl<i V 7 J 

CO M Led: sup sup (-5(e)v^ £ \y,Lo;p) -H(p)) > A5(e)^ |log<5(e)|H . 

r,>0 0<s< V { \y\<CT/e\{e)\p\<L^ ' J 

If we choose A large enough, then (I5.39|) yields that the last event is of probability zero, where we 
have used n = 2 in (jESgP and the fact that e\{e) > 8(e) 2 for small enough e. We deduce that 



liminf inf inf !!^^ > —T 
e^O xeB T 0<t<T X(e) 



1. 



This completes the proof of (i). 
To prove (ii), we define instead 

A(e) := Ae n \loge\*, 

with 5(e) the same as above and A > 1 to be selected, we notice that, for large enough A, 

X(e)>c {A5(ef |log 5(e) 1/{1 ~ 2a) = cA h >(i-25) § ^a | log > M ^a | log § ^ ^ 

where a and (5 are as defined in (12.22)) . We proceed similarly as above, using (|6.7p to obtain 

P) II ueH: sup sup (u e (x,t,w) - u(x,t)) > C 22 TX(e) > 

,>0 0<7<J xeB T 0<t<T J 

CO II Lll: inf inf ( -5(e)v S{e) (y,u ;p) - H(p)) < -X(e)\ 

C p y L G O : inf inf (-5(e)v s ^ (y, u ;p) - H(pj) < -A5(e) a \log5(e)f) 

',>0 0<Y<^ \v\<CT/eX(e)\p\<L\ I j 
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According to (|5.40p . the last set on the right has probability zero provided A > 1 is chosen suffi- 
ciently large. It follows that 

V u e (x,t r ) -u(x,t) 
limsup sup sup — -r < i 

e^O xeB T o<t<T A(ej 

which completes the proof of (ii) . □ 

7. Convergence rates in almost periodic environments 

We conclude by showing that the techniques in the previous section may be used to obtain 
convergence rate for the homogenization of (jl.ip in almost periodic (and in particular periodic) 
media. Since it is necessary to reproduce the qualitative homogenization theory from scratch, we 
take the opportunity to efficiently reorganize and quantify the argument. 

Departing from the hypotheses in the rest of this paper, here we consider H £ C(M. d x R d ) 
satisfying, for each K > 0, the regularity assumption 

H is uniformly continuous on B K x R d and : y G E d | is bounded in C 0,1 (B K ) (7.1) 

and the coercivity condition 

lim inf H(p,y) = +oo. (7-2) 

|p|— >oo y£R d 

The assumption of almost periodicity is that the family of translations of H in the y-variable is 
precompact in the uniform topology of Bk X R d - Precisely, we assume that, for all K > 0, 

{#(•, • + y) : y € R d } is precompact in C{B K x R d ). (7.3) 

We remark that, in this section, we make no convexity assumption on H, nor do we assume any 
analogue of or (i2~Tll . 

The homogenization of coercive Hamilton-Jacobi equations in almost periodic environments was 
proved in [IT]. The key observation was an elegant proof of the following fact. 

Proposition 7.1 (Ishii [UJ). Assume that H G C{R d x R d ) satisfies (f7~I) . fT2D and CL3]). Then 
there exists H £ C(lR d ) such that, for every p £ R d , 

lim sup \5v s (y;p) + H(p)\=0. (7.4) 

The rate of convergence in almost periodic environments follows from Lemma [6 . 1 1 once a rate for 
the limit (|7.4p is obtained. For the latter, it is necessary to quantify the almost periodicity of H, 
which leads us to introduce, for each K > and R > 0, 

Pk(R) '■= sup inf sup \H(p, x + y) — H(p, x + z)\ . 

1 zeB R ( P ,x)eB K xR d 



It is immediate from (|7.ip that, for each K > 0, p/^ is continuous and we see from its definition 
that it is decreasing. The assumption (|7.3p is equivalent to the statement that, for each K > 0, 

lim p K (R) = 0. 

R— >oo 

We next define, for each K > and < 8 < 1, 

77*0*) :=4inf{s>0 : (^) < *} ■ (7.5) 
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The properties of px yield that t\k is a modulus, i.e., for each K > 0, 

t]k '■ (0,1) — > [0, oo) is continuous, increasing and lim7//<-(5) = 0. (7-6) 

5— >o 

Observe that if y i-> H(p,y) is 1-periodic, then Pk{\) = for all K > and, hence, T]k{$) < 5^- 
We also define, for each K > 0, the quantity 



L = L(K) := sup { |g| : inf H(q, y) < sup H(p, y) } , (7.7) 

[ y €R (p,y)£B R xR d J 

which has the property (this is not difficult to check using similar arguments as in Appendix |B|) 
that, for every p G M d , 

\Dv 5 (-;p)\ <L(\p\). (7.8) 

We prove next a quantitative version of Proposition 17. 11 The argument is inspired from |21^ 117]. 
Here we simply reorganize and quantify it. 



Proposition 7.2. Assume that H £ 



H G 



such that, for every K > 0, 5, 7 G (0, 1] and p G -Br-, 



sup 



satisfies (17. lj) . (17. 2|) and (17. 3ft . T/ien iaere existe 



(7.9) 



5«*(y;p) + fr(p) <Vl(S), 



where L = L(K) is given by (|7.7p . 

Proof. The result follows from two facts. The first is that, for every 5 > 0, 

osc5v s (-;p) < r] L (8), 

R d 

and the second is that, for every 5, 7 G (0, 1], 

inf7u 7 (-;p) < sup <!>?/(• ;p). 



(7.10) 



(7.11) 



To demonstrate (|7.10p we fix y G M. d . Let R > be selected below, and choose z G such that 



sup \H{q,y + y) - H(q,y + z)\ < p L {R). 

(q,y)£B L xR d 

It follows from Proposition IB. 1\ by comparing v s (- + y;p) to w' 5 (- + z ;p) ± pl(R)/5, that 



(7.12) 



sup 



6v s (y + y;p) -5v s (y + z;p) < p L (R). 



In particular, and in view of (|7.8p . we have 

sup 6v 5 (y;p) -5v 5 (0;p) < p L (R) + 5L\z\ < p L (R) + 5LR. 



Optimizing over R leads to the choice R := n^((5)/45L, which yields, in light of (|7.5|) . 

tf«*(l/;p)-*«*(0;p)| < \til{S). 

Since y G M. d was arbitrary, we obtain (|7,10p . 

The proof of (|7.1ip is easy and comes from [21]. Suppose, for some 5, 7 G (0,1], that (|7.1ip is 
false. For < a < 1 to be selected, consider the function 



w(y) := v S (y;p) - a (l + |y| 2 ) 
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Using ()7.ip . we see that, if a > is chosen suitable small, then 

H(p + Dw,y) < -sup <$«*(• ;p) + Ca < -inf 7 v 7 (-;p) < H(p + Dv' y ,y) in R d . 

R d Rd 

The comparison principle (Proposition I3.4| ). ()5.2[) and ()T. 1[) imply that, for every R > 0, 

w(0) - w 7 (0) < max(w - t> 7 ) < C f - + - ) - aR. 

Send R — > +co to obtain the desired contradiction. 

It is immediate from (|7.10p and (|7.1ip that, if we define 

H(p) := limM-5v s (Q;p), 

<5->0 

then (I7.9p holds. The continuity of H is immediate from Proposition IB. 91 □ 

The combination of Proposition 17.21 and Lemma 16.11 yields the following convergence rate for the 
homogenization of (jl.ip in almost periodic media. In order to apply Lemma 16.11 we note that its 
proof did not depend in any way on the random environment or the structural assumptions, such 
as level-set convexity or (|2.9p . which are not in force in this section. 



Proposition 7.3. Assume that H G C(R d x R d ) satisfies (fTTT]) . (f73]l and (i7T3j) . Consider the 
unique solutions u £ ,u € C 0,1 (]R a! x [0,T]) of 

uf + H (pu £ , ^) = and u t + H(Du) = in R d x (0, T) 

subject to the initial condition u e (0,t) = u(0,t) = uo(x) 6 C 0,1 (M d ) ; and let K > be such that 

\u £ (x, t) — u s (y, s)\ V \u(x, t) — u(y, s)\ < K(\x — y \ + \t — s\) . 
Then there exists a constant C23 > such that, for all T > 1 and e > 0, 

sup \u e {x,t) -u(x,t)\ < C 23 t(^ +r] L (e^)) , (7.13) 

(x,t)m d x[0,T] v y 



where L = L(K) is given by (|7.7p and i/ie modulus «l(-) fry (|7.5p . 

Proof. For e, a > 0, we define 5(e) := £3 and 

X(e,a) := ((C 22 /T)^(e)) V^)+a 

and observe that (|6.5p holds for a,e > small enough. An application of Lemma 16 . 1 1 yields 

sup \u £ (x, t) — u(x, t)\ < C22TX(e, a). 

(x,t)m d x[0,T] 

Let a -> to get (f7TT3|) . □ 

Observe that for a periodic Hamiltonian satisfying (17. ip and (|7.2p . Proposition 17.31 gives a rate 
of convergence of 0{es) for homogenization. 
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Appendix A. The metric problem: basic properties 

We summarize some elementary facts concerning solution of the eikonal-type equation 

H{Du,y,u>) = n, (A.l) 

in exterior domains, with ji > H*, where H* is a critical parameter. Most of what follows is 
well-known and can be found for example in |20j or [3], but we briefly sketch the arguments for the 
convenience of the reader. Throughout this section, we assume that H satisfies (|2.10p . 

The critical parameter is the infimum of all \x for which (|A.ip admits a global subsolution. 
It turns out (see [3J) that = mm.H and the assumption (I2.9P implies that if* = H(0) = 0. 

We begin with a comparison principle, which makes minimal assumptions on the growth of the 
subsolution and supersolution at infinity. 

Proposition A.l ([4j Proposition 3.1]). Let /j, > 0, K C R d be compact and u, -v € USC(R d \ K) 
satisfy 

H(Du,y,ui) < n < H(Dv,y,uj) inR d \K and limsup (u(y) - v{y)) < (A.2) 

y->K 

and 

liminf#>0. 

Is/I— >oo \y\ 

Then u<v in R d \ K. 

Next we define, for each M. d , fi > 0, and ui E f2, 

m^z, x,ui) := sup {w(z) — w(x) : w € £ and H(Dw, y, u>) < /i in R^} . (A. 3) 

Note that, due to (|2.9p . the zero function belongs to the admissible class, which is therefore 
nonempty. Lemma (|3.5f) and (12. 6p yield that m^y, x, ui) is finite and, in fact, nonnegative and 
bounded from above by C\y — x\, for some C > depending on an upper bound for \i. 

It is immediate from (|A.3P that is measurable with respect to r B®'B®3', since the expression 
on the right of (|A.3P is. Moreover, from (|A.3[) and (|2.2p we see that m M is jointly stationary in its 
first two variables, i.e., for every x,y, z € R d and u) € Q, 

m^(y,x,T z u) = m^{y + z,x + z,u>). (A.4) 

Also immediate from (IA.3P (and the fact that a supremum of a family of viscosity subsolutions 
is a viscosity subsolution, see [HE]) that the m (U (-, x, w)'s are global subsolutions of (lA.ip . i.e., for 
every > 0, i £ R rf and 

H(Dm^(-,x,uj),y,u}) < fj, in R d . (A.5) 

Further properties of the m^'s are recorded in the next proposition. Since detailed proofs of most 
of these facts can be found in [4j, here we present only brief sketches of the arguments. 

Proposition A.2. For every fx > and ui G CI, the following hold: 

(i) For each fixed x € R, the function m At (-,x,w) is a solution of 

H(Dm fl (-,x,uj),y,uj) = fi inM d \{x} and m^(x, x, ui) = 0. (A. 6) 

Moreover, if \x > 0, then is the unique nonnegative solution of (|A.6j) . 
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(ii) If U CI" is open, x € R \ U and u € L is a subsolution of (jA.ip in U , then 

u — mJ-, x, oj) < max (u(y) — mJy, x, oj)) in U. (A. 7) 

y£dU 

(iii) For x,y,z G R d , 

m^(y,x,oj) < m^(y,z,oj) + m^(z,x,oj). (A.8) 

(iv) There exist Z^, > satisfying, for some C, c > depending only on an upper bound for \x, 

ctx<l^<L^< C, (A.9) 

such that 

ln\y - x\< m^(y,x,oj) < L^\y - x\. (A. 10) 

(v) For every x,y £ R rf , 

\m^(y,x,oj) - m fM (z,x,oj)\ < L^y - z\. (A. 11) 

(vi) For every open set U C R d , x € U and y G R d \ {/, 

mJy, x, w) = min (mJy, z, oj) + mjz, x, oj)) . ( A -12) 

z£dU 



(vii) For every x € R rf ; 



H(-Dm^(x,-,uj),y,u) = fi in R d \ {x} (A.13) 



and 

H(-Dm^(x,-,oj),y,oj) < fi inR d . (A. 14) 

(viii) There exists a constant c > 0, depending on an upper bound for [/,, such that, for every 
< // < // and x, y € R d , 

m^(y,x,oj) + c(/i - - y| < m^(y,x,oj). (A. 15) 

Proof, (i) Perron's method yields that m^(-,x, oj) is a solution of (|A.6|) in R d \{x}. Recall from (jA.5P 
that it is also a subsolution in R d . The uniqueness for {i > follows from Proposition lA.il 

(ii) Assume that (|A.7|) fails and, by adding a constant to u, that the right side of (|A.7j) is negative 
while the left side is positive at some point yo € U. Define 



v(y) :-- 



m,j,(y,x,ui) if y<ER d \U, 

u{y)y m^(y,x,oj) if y € 17, 



and observe that i>(x) = 0, m^,(-,x,oj) < v in R d and m^(yo, x,oj) < v(yo). Moreover, it is clear 
from its definition that u is a subsolution of (jA.ip in R d . This contradicts (|A.3p . 

(iii) follows from (ii). Indeed, observe that, according to Lemma l3,5j as a function of y, the left 
side of (|A.8j) is a subsolution of (jA.ip in R rf and is equal to the right side of (I3.5P at y = z. The 
maximality of m At (-, z,oj) (i.e., property (ii) with U = R d \ {x}) then implies (IA.8p . 

(iv) The lower bound of (jA.lOp follows from the observation that, if \i < K, then (|2.5[) and (j2.9|) 
imply the existence of some c > 0, depending on K, such that, for every x € R rf , the function 
y i — y cfi\y — x\ is a subsolution of (jA.ip in R rf . The upper bound is immediate from Proposition (jA.ip 
and the fact that, for large enough C > and any x € R d , the map y i-> C|y — xj is a supersolution 
in R d \ {x}. 

(v) The Lipschitz estimate (jA.lip is immediate from (jA.8h and (jA.lOj) . 
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(vi) One direction of (|A. 12() is obvious from (|A.8|) and holds without restriction on x, y G R d and 
uj G 0. For the other, we first assume that U is bounded and note that, due to Lemma 13.51 

(j>(y) := min (mJy, z, u) + mJz, x, u)) 

zGdU 

is a solution of H(Dcj),y,oj) < fj, in R \ f7. For y G 9C7 we may take z = y in the minimum to 
obtain 4>(y) < m^(y, x,lo). The maximality of m^{-,x,uj) yields <fi < mA-,x,oj) in R d \ U, which 
is the other side of (lA.12p . If U is not bounded, then we approximate U by bounded sets and 
use (jA.lOp . which implies that points z G dU which are very far away from x and y are irrelevant. 

(vii) follows from Lemma 13.51 and (|A.3p , Indeed, it is immediate from the definitions that 
m^(y, x,u) = ria(x,y,u}), where is the function corresponding to for the Hamiltonian 
H(—p,y,uj). Thus (|A.13P holds and we can apply Lemma 1331 to obtain (|A.14jh 

(viii) follows from (|2.5p and (ii). We observe that, as a function of y, the left side of (|A.15j) is 
a subsolution of (jA.ip . provided that we take c > small enough. Both sides of (|A.15|) vanish at 
y = x, and so the result follows from the maximality of m»(',i, cj). □ 

Remark A. 3. The functions can be expressed by the following representation formula due to 
Lions [20], which provides the above facts with a control theoretic interpretation: 

mJy,x,u)) = M I JJ r Y l (s),i(s),uj)ds : 7 eC(x,y)} , (A.16) 



where Q(x,y) is the set of Lipschitz curves 7 : [0, 1] -> R d such that 7(0) = x and 7(1) = y and J, 



f 



is the support function of the /i-sublevel set of H, given by 

J/j,(.Q, V, := sup {p • q : H(p, y, u) < fj} . 

The expression (lA.16j) provides us with an interpretation of m^(y, x,u) as measuring the "cost" 
of moving from the point x to the point y in the medium u. We make no direct use of ()A. 161) in 
this paper, preferring instead to work with the maximality property (Proposition IA.2l fii)) which 
is equivalent to it. Nevertheless, our intuition is enriched from (|A.16P and it suggests an analogy 
between the metric problem and first-passage percolation. 

Definition A. 4. For each [i, t > and x G R rf , we define the reachable set to x in time t by 

^t(x) := {(j/.w) G R d x O : m„(y,x,oj) < t} . (A.17) 

We also use the notation $.^ t = 3l M)t (0), 0^ t {x) := {y G M d : (y,w) G Ol^ t (x)} and Dl£ t = ££ >t (0). 

We continue by examining some elementary properties of the reachable set. According to (jA.lOp . 
we have m^(y, x,u)) < t for y in the interior of Jl^ t (x) and, for all fj,, t > and u; G fJ, 

%^^ t CB (/v (A.18) 

In fact, from (lA.lOh and (lA.20p . we have the following improvement of (IA.18p . which estimates the 
rate of growth of the reachable set: for every < s < t and w£fl, 

{x G R d : dist < L'H* " *)} Q X^,t C {x G R d : dist (x.tt^) ^ " s )} • ( A - 19 ) 

It is useful to note that, in the particular case that x = and J7 = 3?^, Proposition IA.2( vi) asserts 
that, for every t > and y 

mJy,0,u)) = t+ min mJy,z,u). (A.20) 
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It also follows from (|A,12p that, for every s, t > and y G M. d , 

m^(y,^ t ,w)-m M (i/,^ SJ w) < |s-t|. (A.21) 

The maximality property (Proposition I A. 2Y ii)) can be improved for the domain U = %!£ t (x)\{x} 
by restricting the maximum over dU to {x}, as stated in the following lemma. This is the crucial 
fact that localizes the metric problem. 

Lemma A. 5. For every oj G ft, \i > 0, x G M. d and w G Lip(3?^ )t (x)), 

H(Dw,y,uj) < fi in 0l^ t (x) implies w(-) — w(x) < m^(-,x,oj) in (A. 22) 

Proof. Let ui G Lip(3^ 4 (aj)) satisfy the antecedent of (jA.22j) and assume with no loss of generality 
that w(x) = 0. Define w := m^(-,x,oj) V (w(-) A (t — e)) and observe that, in light of Lemma [331 
and (|2.9p . u; is a subsolution of H(Dw,y,oj) < fi in We deduce from the maximality property 
that w < m^(-,x,oj) in R rf and in particular w A (t — e) < m /1 (-,x,a;) in 3?^ t (x). Sending e — > 
and using that m^(y,x,oj) < t in the interior of 3^ t (x), we obtain that w < m^(-,x,oj) in Jl^ t (x), 
as claimed. □ 

It follows from Lemma IA.5I that the representation formula (|A.3|) may be restricted to the 
reachable set. Precisely, for every oj G fl, > 0, x G M rf and y G 3J^ jt (x), 

m IJi (y, x, oj) = sup{if(y) — u>(x) : w G Lip (3?" t (a;)) and H(Dw,y,uj) < \i in 3?^ ,t(^)} • (A. 23) 

This is immediate from (|A.22|) . 

In the proof of Lemma |3.1| we require a refinement of Lemma lA.51 To this state this, we define, 
for every nonempty closed set K C M, d , x, y G K and oj G ft, 

m„ (y, x, oj) := sup{u>(y) — w(x) : io G Lip (if) and H(Dw,y,oj) < \i in if}. (A. 24) 

It is immediate that m^f is S(-?0-measurable, and 

K\ C if 2 implies that m^ 1 (y,x,oj) > m^ 2 (y,x,oj) for all x,y E K\. (A. 25) 

Thus (|AT23|) yields that 

R^ t ^KQR d implies that m^; 0, w) = (•, 0, w) in ft^. (A. 26) 

In fact, by a similar argument as in the proof of Lemma lA.51 if G K, then 

inf m^(y,0,oj) > t implies that m5(-, 0, oj) = mJ-, 0, oj) in {m^(-, 0, oj) < t\ . (A.27) 

y£dK 

Here is the proof: consider the function mff (•, 0, oj) A (t — e) and extend this to M. d by giving it the 
value (t — e) outside of K. This is a global subsolution by Lemma 13.51 and the assumption that 
m^(-,0,oj) >tondK. Thus by the maximality of we deduce > mjf(-, 0, oj) A (t — e). Sending 
£ — y yields m At (-,0,o;) > m^(-,0,oj) in |m^(-,0,a;) < i}. In light of (|A.25|) . this completes the 
proof of (|A.27p . which is needed in the proof of (|3.7|) . 
For /j, > 0, we define the filtration {S'^, i t}t>o by 

3^ fi = {0,fl} and S^t := S(3W) for t > 0, 

where 9 is as defined in Section [2j Intuitively, the cr- field 3ut consists of the information uncovered 
by observing the Hamiltonian H(p,x,oj) for x G %tt(u;) and p G M d . 
It is immediate from (|A.23p that, for every y G M. d , 

oj mfj,(y,0, u})t^ ye3tiitt y{u)) is ^-measurable. (A.28) 
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According to (|A.18|) and (|A.28|) . for every y G R d and t> L^\y\, 

uj i—)- m^(y,0,uj) is 3~^-measurable. (A. 29) 

We next define, for every compact set K of M. d , y G M. d and uj G f2, 

mJy,K,u)) := inf mJy,z,uj) and mJK,y,uj) := inf mJz,y,u>). (A. 30) 

Proposition A. 6. T/ze functions defined in (|A.30p are f/ie unique nonnegative solutions of 

H(Dm^(-, K,uj),y,u) = \i and H(—Dm tl (K,-,u),y,uj)=fJL inR d \K (A. 31) 

which satisfy the boundary condition m^(-,K,uj) = m^(K, ■, uj) = on dK . Moreover, they are 
also given by the representation formulas 

mJy, K, uj) = min sup I w(z) — w(x) : w G £, H(Dw, z,w) < u inM d \K> , (A. 32) 

and 

mJK, z, uj) = min sup \ w(x) — w(z) : w G L, H(Dw, y,oj) < fi in M. d \ K > . (A. 33) 

Proof. That m jU (- ; X, cj) and m^(K, ■, uj) satisfy ()A.3ip is immediate from (|A.6P and ()A.13P and the 
fact that, in light of Lemma 13.51 the infimum of a family of solutions is a solution. Uniqueness 
follows from Proposition lA.il 

It is easy to check from the standard Perron argument, using again the fact that a minimum 
of solutions is a solution, that the right side of (|A.32j) . as a function of z, is a solution of the 
first equation in (IA.31j) . Clearly it vanishes on K by taking x = z in the minimum and the zero 
function in the supremum. Therefore it must equal m^(-,K, uj) by the uniqueness of the latter. 
This proves (IA.32D . and (|A.33h is obtained similarly. □ 

It is immediate from (|A.32p and (|A.33P that, for every compact K 6 R rf and y G W 1 , 

uj i— > m^(y, K,uj) and w H m (1 (A',j/,w) are S (M d \.fC)-measurable. (A. 34) 

The connection between the metric problem and H is provided by the equivalence 

fi>H(p) if and only if liminf (m^y, 0, uj) — p ■ y) > a.s. in uj. (A. 35) 

|y|— >oo 

This fact is at the heart of the qualitative homogenization theory. See [4j for a proof. Using (jA.8j) 
and the subadditive ergodic theorem, we can find a positively homogeneous and convex function 
rrifi G C(M rf ) such that, for every y,z G M. d and uj belonging to some fixed event of full probability, 

lim - \mJty,tz,uj) - mJy - z)\ = 0. (A.36) 

t— >oo t 

For the details of the argument, we refer to [3j 0]. It follows from ()A.8P and (|A.36|) that m M is 
convex and positively homogeneous. From ()A.35P we deduce that (c.f. [1]) 

H(Dm^) = ii in R d \ {0}. (A.37) 

The functions prescribe the effective Hamiltonian H{p) via the formula 

H(p) := inf {/i > : m M (y) > p • y for all y G M d } , (A.38) 

a fact which follows easily from ()A.35P and ()A.36P and is essentially equivalent to (|A.37p . We 
deduce from (|A.37|) that H is level-set convex. Moreover, it is a consequence of (12. 7p and (12. 8p 
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that, for every fi > 0, the level set {p G M. d : H(p) = /i} has empty interior. Of course, the minimal 
level set of H may indeed have nonempty interior. 

Finally, from elementary convex geometric considerations (see [3] once again) we deduce that, 
for every p int{g : H(q) = 0}, we can find a direction e so that the plane with slope p touches 
from below at w. Precisely, setting /x := H(p), there exists e G M d with |e| = 1 such that 

mJe) — p ■ e = = min (mJy) — p-y) . (A. 39) 

y£R d 

Appendix B. The approximate cell problem: basic properties 

We summarize the properties of the approximate cell problem 

5v s + H(p + Dv 5 ,y,Lo) = in R d . (B.l) 

Here p G W 1 and 5 > are given parameters. We note that the results in this section do not make 
use of assumptions (|2.8p or of (|2.9p . nor do they depend in any way on the random parameter u, and 
therefore they hold (with appropriate changes in the notation) for the Hamiltonians encountered 
in Section 

We begin with a comparison principle for (|B.ip . which can be reduced to Proposition 13.41 by an 
argument which perturbs the solutions by adding appropriate terms with linear growth (alterna- 
tively, a proof can be found in [HE]). 

Proposition B.l. Let 5 > 0, p G ~R d and wed. Suppose u, —v € USC(M d ) satisfy 

5u + H(p + Du,y,ui) < < 5v + H(p + Dv,y,u) inR d , (B.2) 
and v is bounded below on M. d . Then u < v inM. d . 

For each 5 > 0, p G M. d and w € O, we define 

v s (y, oj;p):= sup lw(y) : w G USC(M d ) satisfies 5w + H(p + Dw, y,oj)<0 in M d | , (B.3) 

with the differential inequality interpreted either in the viscosity or in the almost everywhere sense, 
as these are equivalent by Lemma 13.51 It is clear that the constant function 

w := — - ess sup H (p, y, oj) = — - ess sup H (p, 0, •) 

belongs to the admissible class, hence v s (-,oo ;p) > w. Similarly, the function 

v := — — essinf H(p, 0, •) 

is a bounded supersolution of (jB.ip . and Proposition IB. II yields v s (-,u;p) < v. We have shown 
that, for all y,p G M d , ui G O and 5 > 0, 

— esssup H(p, 0, •) < 5v s (y,co ;p) < — essinf H(p, 0, •). (B.4) 
n n 

Immediate from (|B.3P and (|2.2p is that the v s, s are stationary functions. That is, for all y,z,p G M. d , 
and 5 > 0, 

v 5 (y,T z u;p) = v S (y + z,io;p) (B.5) 
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We summarize some further properties of the v s, s in the following proposition. The proofs, which 
are standard in the theory of viscosity solutions, are only sketched. For the statements we need to 
define the constants 



K p : = sup < \p — q\ : ess inf H (q, 0, •) < ess sup H (p, 0, •) > (B-6) 
{ n n J 

and 

lip := max sup ess sup \H(p + q, 0, •) — H(p + (1 + a)q, 0, •) I . (B.7) 
ff=±1 kl<^p ^ 

Note that K p is finite for \p\ bounded by (|2.6p and U p is also finite for bounded \p\ by ()2.5|) . 
Proposition B.2. For every p € M. d , U) G $7 and 5 > 0, the following hold: 

(i) The function v s (•, uj ; p) is the unique solution of (jB.ip belonging to BUC(lR d ). 

(ii) For every x,y € M. d , 

\v 5 (x,Lo;p) -v 5 (y,uj;p)\ < K p \x - y\. (B.8) 

(iii) For all q G R d , 

\5v s (y, uj ;p) — 5v s (y, oj ; q)\ < sup ess sup \H(p + z, 0, •) — H (q + z, 0, •)! . (B.9) 

|z|<ftrpV^ g n 

(iv) For every n > 8, 

\Sv\y,uj; P )-riv*(y,u>; P )\<n p (l-?) . (B.10) 

Proof, (i) The boundedness of v s was proved above in (1B~4|) . That (IB~3]l gives a solution of (IB.lj) in 
BUC(M d ) follows from the classical Perron argument adapted to viscosity solutions (see [9]). The 
uniqueness of v s is immediate from Proposition lB.il 

(ii) According to (|B.4|) and Lemma 13.51 \p) is Lipschitz with a uniform constant that can 
be estimated by (|2.6p . We end up with (|B.8p with K p as defined in ()B.6P . 

(iii) The dependence of 5v s on p can be controlled using the comparison principle together 
with (|2.5|) and (|B.8p . The argument is routine, so we merely sketch it. One inserts v s (-,oj ; q) into 
(jB.lj) . adds or subtracts a constant until the resulting function is a supersolution or subsolution, 
and applies Propostion lB.il The estimate produced by this argument is (|B.8[) . We remark that, due 
to (12. 5p and (12. 6|) . the right side of (|B.9p is controlled by C\p — q\ for a constant C > depending 
on an upper bound for \p\ V |g|. 

(iv) The dependence of 5v s on 5 is also controlled with a simple comparison argument and the 
help of (|2.5p and (]B.8p . Set A := 5/r/ and w(y) := Xv 6 (y,u> ;p) and check that 

nw + H(p + Dw, y, u) < Sv s + H(p + XDv s , y, oj) < C(l - A) in R d , 

where C = H p defined in (|B.7p . An application of Proposition IB. II yields 

5v 5 (;u);p) = r]w< nv* 1 (■ , uj ; p) - C(l - A), (B.ll) 

which is half of (IB.lOp . The other inequality is obtained via a similar argument. □ 

We conclude by recalling that the connection between (IB.ip and H lies in the limit 

suplimsup sup \Sv s (y, u ;p) + H(p)\ = 0, (B.12) 

R>0 <5-s>0 yeB R/s 
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which holds for all p £ M. d and uj € Qo, where Oo £ 3~ is a event of full probability which does not 
depend on p. This assertion is equivalent to homogenization (see [5j for a proof). 
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